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INTRODUCTION 


The  SDCS  attenuation  experiment  was  motivated  by  indications  of  differ¬ 
ences  in  the  m^-Ms  relationship  for  explosions  at  NTS  relative  to  explosions 
at  USSR  test  sites.  These  differences  can  be  explained  by  anomalously  high 
attenuation  in  the  mantle  under  the  various  test  sites  in  the  Western  United 
States  as  opposed  to  a  high  Q  mantle  under  shield  areas  of  the  USSR.  This 
report  consists  of  three  parts  and  summarizes  all  of  the  results  of  the 
project  to  date. 

The  first  part  discusses  the  analysis  of  data  from  the  SDCS  station 
network.  This  data  base  is  also  supplemented  by  the  analysis  of  some  LRSM 
data  in  order  to  extend  the  areal  coverage  and  to  compensate  for  the  fact 
that  only  one  SDCS  station  was  located  on  a  shield. 

The  second  part  consists  of  a  detailed  discussion  of  the  results  of  the 
SDCS  experiment  in  the  framework  of  research  conducted  on  attenuation  by  us 
and  others.  This  section  outlines  regional  variations  of  Qa  and  Qg  under 
the  United  States  and  puts  limits  on  the  absolute  and  relative  variations  of 
these  quantities.  The  need  for  a  frequency  dependence  of  Q  is  also  dis¬ 
cussed.  Since  no  comparable  study  of  short-period  attenuation  has  been 
undertaken  in  the  past  and  since  short-period  data  are  influenced  by  many 
extraneous  factors  not  related  to  mantle  Q,  a  methodology  had  to  be  created 
to  reduce  the  effects  of  these  extraneous  factors  and  to  choose  parameters 
less  sensitive  to  them  and  more  diagnostic  of  mantle  Q.  The  data  are  inter¬ 
preted  in  the  following  framework: 

a)  Due  to  the  well-documented  focusing  effects  of  small  scale 
inhomogeneities  in  the  crust  and  the  uncertainties  in  the  estima¬ 
tion  of  crustal  amplification,  short-period  body  wave  amplitude 
anomalies  are  likely  to  be  biased  and  are  not  reliable  indicators 
of  mantle  Q. 

b)  While  short-period  body  wave  spectral  shapes  (especially  falloff 
rates  at  high  frequencies)  are  extremely  sensitive  to  small  varia¬ 
tions  of  mantle  Q,  they  are  considerably  less  affected  by  other 
factors. 
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c)  To  measure  mantle  attenuation  a  combined  interpretation  of  P  and  S 
wave  amplitude  and  spectral  data  is  necessary,  with  more  weight 
being  given  to  spectral  measurements. 


The  third  part  of  this  report  contains  our  arguments  for  using  the 
methods  outlined  above.  We  feel  that  justifying  the  methodology  within  the 
main  body  of  the  report  would  continually  divert  the  reader's  attention  to 
side  issues;  therefore,  we  have  reserved  a  separate  part  of  the  report  to 
demonstrate  the  validity  of  points  a-c  above  and  to  discuss  several  other 
factors  that  need  to  be  addressed  to  overcome  possible  objections.  This 
is  an  essential  part  of  the  report  and  could  be  read  first.  Part  III 
consists  of  Sections  A-F,  each  of  which  is  referred  to  in  the  main  part 
of  the  text. 


In  this  report  we  describe  attenuation  in  terms  of  the  quantity  t*  =  J  — , 
where  Q  is  the  quality  factor  along  a  certain  seismic  ray  path  and  t  is  the 
travel  time,  t*  is  naturally  a  path-dependent  quantity  since  it  is  a  function 
of  the  Q  variation  in  the  Earth  which  in  turn  is  a  function  of  depth  and 
region.  It  is  probably  also  frequency  dependent,  t*  appears  to  be  the  most 
convenient  parameter  for  characterizing  attenuation  because  for  distances 
greater  than  25  degrees  it  changes  little  with  epicentral  distance.  Thus,  it 
can  be  used  as  a  single  (although  frequency  dependent)  parameter  to  charac¬ 
terize  regional  differences  for  various  types  of  paths  described  in  terms 
of  upper  mantle  structures  under  the  source  and  receiver.  For  distances  less 
than  20  to  25  degrees,  t*  increases  with  epicentral  distance  due  to  the  fact 
that  the  body  waves  progressively  penetrate  the  low-velocity  -  low-Q  layer 
(LVZ)  in  the  upper  mantle  as  the  distance  increases.  The  contributions  from 
the  LVZ  (which  is  coincident  with  the  low  Q  layer)  dominate  the  integral  and 
thus  t*  becomes  a  property  that  depends  almost  entirely  on  the  types  of  upper 
mantle  structures  the  ray  path  crosses.  Thus: 


“■/V  /V?  (1) 

D  ^  U  M 

where  D  is  the  integral  along  the  downgoing  leg  of  the  ray  path  through  the 
LVZ  and  U  is  the  upgoing  part,  t*  is  the  contribution  from  the  rest  of  the 
path  through  the  crust  and  the  lower  mantle.  It  is  relatively  small  because 
the  crust  and  lower  mantle  are  high-Q  regions. 
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Alternate  representations  of  attenuation  are  leas  convenient.  For 
instance,  if  the  apparent  Q  is  used  it  has  to  be  specified  as  a  function  of 
travel  time  and  other  parameters.  It  should  be  understood  therefore  that 
t*  in  this  report  always  refers  to  values  measured  at  teleseismic  distances 
unless  it  is  clearly  associated  with  near  distances. 

Throughout  this  report  the  mean  amplitude,  travel  time,  t*  and  dominant 
period  differences  between  stations  are  not  considered  significant  unless 
they  exceed  the  95%  statistical  confidence  levels.  The  confidence  level  of 
some  of  the  data  is  higher  than  99%.  This  kind  of  stability  is  not  often 
found  in  geophysical  literature,  and  the  existence  of  interstation  differen¬ 
tials  in  the  quantities  measured  can  hardly  be  questioned  in  such  cases. 

Details  of  the  experiment  that  are  not  essential  to  the  main  conclusions 
of  the  report  are  presented  in  appendices  at  the  end. 
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PART  1 

DATA  ANALYSIS  FROM  SDCS  AND  SELECTED  LRSM  STATIONS 

The  basic  theory  of  the  SDCS  experiment  is  that  by  measuring  relative 
amplitudes  and  spectral  ratios  for  common  events  between  stations  one  can 
infer  the  degree  of  anelastic  attenuation  under  each  station.  The  design  of 
the  experiment  by  DARPA  is  based  on  an  implicit  assumption  of  a  degree 
of  seismic  reciprocity. 

Seismic  reciprocity  is  a  general  law  stated  by  Knopoff  and  Gang!  (1959), 
Hudson  (1969),  and  others  and  is  valid  for  arbitrary  inhomogeneous  linearly 
elastic  media.  This  law  expresses  the  interchangeability  of  source  force 
systems  and  combinations  of  receiver  displacement  parameters.  For  example, 
one  can  state  the  interchangeability  of  dilatational  force  and  displacement 
systems  (Hudson,  1969;  Knopoff,  1979).  Reciprocity  requires  that  the  sources 
and  receivers  when  interchanged  occupy  exactly  the  same  positions  in  space. 
This  condition  is  not  satisfied  in  a  strict  sense  for  the  SDCS  experiment 
at  the  NTS  sites.  This  is  because,  for  the  m^  -yield  application,  the 
sources  are  at  the  explosion  sites  and  the  receivers  are  mostly  WWSSN  and 
LRSM  stations  scattered  world-wide,  whereas  for  the  reciprocal  problem, 
the  receivers  are  close  to  the  explosion  site  (but  not  at  the  same  depth) 
and  the  sources  are  located  at  the  seismic  belts  of  the  earth.  Thus  from 
basic  physical  grounds  alone  one  cannot  expect  reciprocity  to  be  exactlv 
valid  for  the  NTS  experiment. 

The  Special  Data  Collection  System  network  consisted  of  a  variable 
number  of  portable  seismic  stations  deployed  across  the  United  States  and 
Canada.  It  is  essentially  an  extension  of  the  LRSM  (Long  Range  Seismic 
Measurement)  network  with  more  up-to-date  digital  recording  systems,  although 
at  some  sites  analog  recording  was  still  used.  Between  1976  and  1979,  these 
portable  stations  were  deployed  at  various  nuclear  test  sites  in  the  western 
United  States  and  at  the  sites  RKON  (Red  Lake,  Ontario),  IFME  (Island  Falls, 
Maine),  and  HNME  (Houlton,  Maine)  in  order  to  measure  the  differences  in 
anelastic  attenuation  in  the  mantle  under  each  site.  The  sites  are  marked  on 
the  map  in  Figure  1,  and  relevant  site  data  are  tabulated  in  Table  I. 

The  instrument  responses  of  the  short-period  systems  were  identical  at 
all  sites  except  HNME.  This  makes  it  impossible  to  directly  compare  dominant 
signal  periods  and  trace  amplitudes  at  this  station  with  the  rest  of  the 
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Figure  1  Locations  of  the  SDCS  and  LRSM  stations  analyzed  in 
detail  in  Part  I  of  this  report. 


Description  of  SDCS  and  LRSM  stations  discussed  in  Part  I  of  this  report 
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SUBSCRIPTS  INDICATE  ORIENTATION  OF  R,  T  COMPONENTS;  NTS  INDICATES  ORIENTATION  TOWARD  THE  NEVADA  TEST  SITE 


stations.  The  instrument  responses  of  several  kinds  of  systems  discussed 
in  this  report  are  shown  in  Figure  2. 

The  approach  taken  for  event  selection  was  to  measure  events  on  the 
first  published  lists,  which  were  the  bulletins  of  the  Hagfors  Observatory, 
that  of  the  French  national  network  or  the  PDE  lists  of  the  USGS.  Locations 
were  double  checked  with  later,  more  accurate  event  lists.  These  preliminary 
locations  were  sufficient  for  the  purpose  of  relative  magnitude  measurements. 

Results  of  Amplitude  Studies 

During  the  course  of  the  SDCS  experiment  the  amplitudes  and  periods  of 
more  than  600  events  were  measured  for  the  purpose  of  computing  interstation 
m^  differentials.  The  events  were  all  in  the  range  25®  <  A  <  85°  from  the 
stations  to  ensure  traversal  of  the  upper  mantle  by  most  of  the  P  wavetrain. 

A  complete  listing  of  the  events  is  shown  in  Appendix  A. 

Since  the  SDCS  experiment  was  designed  to  resolve  problems  arising  from 
a  study  of  Mg-m^,  it  was  necessary  to  compute  the  conventional  using  the 
maximum  amplitude  in  the  first  3  seconds  of  each  signal.  The  writers  are  of 
the  opinion,  in  agreement  with  Butler  (1979),  that  the  amplitude  of  the  first 
cycle  or  "b"  phase  is  more  meaningful;  nevertheless,  we  shall  show  that  our 
results  agree  extremely  well  with  those  obtained  from  "b"  phase  measurements. 
Following  Cleary  (1967)  and  Butler  (1979)  we  have  also  compiled  differentials 
of  trace  amplitudes  corrected  for  distance,  and  we  have  abandoned  the  compu¬ 
tation  of  the  uncorrected  magnitude  m'  used  in  previous  reports  because  we 

a 

feel  that  it  is  less  meaningful  than  the  trace  amplitudes  (see  section  B  of 
Part  3). 

Histograms  of  the  and  trace  differentials  are  shown  in  Figures  3  and 
4  for  the  station  pair  RKON  and  0B2NV,  two  stations  of  great  diagnostic 
value  since  both  of  these  are  located  on  granite  and,  presumably,  are  free 
of  large  crustal  amplification  effects.  These  figures  show  the  extremely 
large  scatter  in  the  amplitude  and  magnitude  differentials  characteristic  of 
short-period  data.  This  is  not  surprising  in  view  of  the  comments  in 
section  D  of  Part  3  of  this  report.  A  noteworthy  data  point,  the  Novaya 
Zemlya  nuclear  explosion  of  02  OCT  76  (marked  as  N.Z.),  is  on  the  flanks  of 
the  distributions  in  spite  of  the  fact  that  the  source  is  axisyraraetrlc. 
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Figure  3  Histogram  of  m,  differentials  between  RKON  and  OB2NV. 

The  mean  m^  differential  is  0.173  ±  .066  (95%  confi¬ 
dence)  magnitude  units.  The  point  N.Z.  marks  a  data 
point  from  Novaya  Zemlya. 
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Ado910  Au) 


Figure  4  Histogram  of  P  wave  trace  amplitude  differentials 

between  RKDN  and  0B2NV  using  ten  base  logarithms  of 
amplitudes  as  units.  The  differential  is  0.267  ± 
.064  (95%  confidence)  showing  that  the  RK0N/0B2NV 
amplitude  ratio  is  about  1.85. 


Histograms  of  the  magnitude  and  trace  amplitude  differentials  for  the  rest 
the  station  pairs  are  collected  in  Appendices  B  and  C  of  this  report. 


Consider  now  the  standard  deviation  of  relative  magnitude  differ¬ 

entials  as  a  function  of  the  distance  between  station  pairs.  Figure  5 
is  a  plot  of  o.  versus  log  A°.  A  linear  empirical  relationship  appears 
to  exist.  The  causes  of  the  linear  increase  in  a.  are:  1)  decreasing 

ATP 

signal  similarity  for  the  P  waveforms  as  the  distance  increases  so  that  m^ 

cannot  be  consistently  picked  at  similar  portions  of  the  P  wavetrain,  2) 

source  radiation  patterns  that  cause  systematic  variations  of  amplitudes, 

and  3)  increasing  effects  of  focusing  and  defocusing.  In  order  to  take 

into  account  the  effect  of  signal  coherence  evidenced  by  the  a,  -A° 

Am, 

relationship,  we  make  use  of  a  least  squares  technique  which  also  allows  us 
to  incorporate  all  of  our  measurements  into  the  determination  of  magnitude 
differentials  across  the  SDCS  network.  For  a  given  event  we  write  the 
magnitude  differentials  between  pairs  of  stations  i  and  j  in  the  form 

AmJ-AmJ«kA^  +eij(^j)  (1) 


i  k  ii 

where  Am^  is  the  station  term  (bias)  of  station  i,  Am^"1  *s  the  observed 

magnitude  differential  between  station  i  and  j  for  event  k  and  is 

an  error  term  dependent  upon  the  distance  A0^  between  the  oca'. ions.  The 

expected  value  of  this  error  term 


B‘‘i,  «*»  ■  w5j> 


(2) 


can  be  read  from  the  regression  line  in  Figure  5  as  a  function  of  distance. 


Because  taking  differences  in  all  possible  combinations  is  redundant, 
a  hierarchy  of  stations  was  assigned  as  follows:  0B2NV,  0B3NV,  YFNV, 

YF4NV,  YF2NV,  YF3NV,  NTNV,  NT2NV,  FANV,  GBNM,  RKON,  HNME,  and  we  used  the 
leftmost  available  station  bias  term  as  the  reference  (positive)  term  in 
equation  (1).  This  hierarchy  was  selected  to  optimize  the  distance 
distribution  so  that  most  of  the  distances  used,  and  consequently  the 
values  of  a.  ,  are  as  small  as  possible.  Thus  0B2NV,  at  the  center  of  the 
NTS  cluster  and  located  between  GBNM  and  FANV,  is  the  prime  candidate  for 
the  most  commonly  used  reference  station,  while  RKON  and  HNME  are  outliers 
and  rank  low. 


This  least  squares  procedure  causes  only  minor  changes  in  the  relative 
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Figure  5  Standard  deviation  of  relative  differentials  vs. 

interstation  distance.  The  increase  reflects  increasing 
variability  of  nt  readings  for  more  distant  station 


®h  and  trace  amplitude  differentials  compared  to  those  directly  computed 
from  data  and  given  in  the  histograms  in  Appendices  B  and  C.  The  means  of 
the  differentials  relative  to  0B2NV  of  and  trace  amplitudes  are  listed 
in  Table  II  and  also  shown  with  their  95%  confidence  limits  in  Figures  6 
and  7.  The  most  prominent  features  of  these  figures  are  the  high  amplitudes 
of  m^  values  at  the  stations  on  Yucca  Flats  (YF's)  and  on  Pahute  Mesa  (NT's). 
The  low  amplitude  of  the  station  IFME  in  Maine  is  another  feature  of 
Figure  6. 

In  order  to  correct  for  the  effect  of  the  crust  on  these  m^  and  trace 

amplitudes  we  compared  changes  in  synthetic  P  pulse  amplitudes.  Pulse  shapes 

for  a  50  kt  explosion  as  modeled  by  von  Seggern  and  Blandford  (1972)  were 

computed  after  passing  them  through  layered  halfspace  models  using  the 

standard  Haskell  algorithm.  Table  III  lists  the  model  parameters  used  in 

the  calculations.  We  attenuated  each  pulse  with  a  multiplicative  spectral 

factor  exp(-irft*),  where  t*  was  chosen  to  be  0.45,  a  typical  value  for  the 

WUS.  By  removing  most  of  the  high  frequencies,  the  attenuation  factor  makes 

the  pulse  more  rounded.  This  pulse  has  a  spectrum  which,  in  spectral  content, 

well  represents  the  average  teleseismic  P-wave  arrivals.  It  is  peaked  at 

-2 

1  Hz  and  falls  off  at  a  rate  of  somewhat  more  than  w  at  high  frequencies. 

To  reduce  variations  caused  by  changes  in  the  angle  of  incidence,  we  com¬ 
puted  the  synthetics  for  three  angles  (20°,  25°,  and  30°  measured  from  the 
vertical)  and  then  averaged  the  relative  amplification  factors  between 
stations  obtained  for  these  three  angles. 

The  Yucca  Flats  sites  rest  on  thick  unconsolidated  sediments  and  tuff 
that  cause  considerable  signal  amplification  (Houser,  1968;  Fernald  et  al, 
1968;  Healy,  1968;  Ramspott  and  Howard,  1975;  Hays  and  Murphy,  1971).  The 
FAULTLESS  site  (FANV)  is  also  located  over  alluvium  and  tuff,  but  according 
to  test  site  Information,  the  alluvium  is  more  consolidated  at  this  site 
(McKeown  and  Dickey,  1969)  (Lt.  Col.  Bulin  of  ARPA  also  provided  us  with 
data  relevant  to  the  FANV  site).  Alluvium  and  a  thick  sedimentary  carbona¬ 
ceous-shale  sequence  also  underlay  the  GASBUGGY  site  (GBNM)  (Thornbrough, 
1971).  In  all  of  these  cases  the  part  of  the  structure  that  primarily 
determines  crustal  amplification  was  found  to  be  near  the  surface.  Conse¬ 
quently,  the  structures  were  modeled  only  down  to  the  basement  (Der  et  al. 
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TABLE  III 


Crustal  models  used  lor  estimating  crustal  amplifications 
at  the  SDCS  stations 


(i(tcm) 

a(km/sec) 

SOcm/sec) 

p(g/cm2) 

FA- MV 

0.44 

2.50 

1.06 

2.30 

1.00 

3.00 

1.60 

2.30 

1.00 

3.50 

1.73 

2.30 

1.525 

4.00 

2.20 

2.50 

as 

5.70 

3.29 

2.70 

GB-NM 

0.700 

2.00 

1.06 

2.00 

0.1S0 

3.25 

1.63 

2.00 

0.200 

3.45 

1.73 

2.00 

0.120 

3.90 

2.16 

2.28 

0.600 

4.40 

2.48 

2.49 

0.625 

4.80 

2.73 

2.55 

00 

5.70 

3.29 

2.70 

CQ-NV 

2.0 

6.10 

3.60 

2.70 

ao 

6. 10 

3.60 

2.70 

HN-ME 

10.0 

5.90 

3.36 

2.70 

“ 

6.35 

3.62 

2.72 

1F-ME 

2.0 

6.10 

3.60 

2.70 

00 

6.10 

3.60 

2.70 

NT-NV 

1.0 

3.00 

1.80 

2.00 

4.0 

3.60 

2.00 

2.20 

5.0 

5.70 

3.36 

2.70 

OD 

6.10 

3.60 

2.80 

NT2NV 

1.0 

2.86 

1.75 

2.00 

4.0 

3.60 

2.00 

2.20 

5.0 

5.70 

3.36 

2.70 

*> 

6.10 

3.60 

2.80 

0B2NV 

10.0 

5.70 

3.36 

2.70 

00 

6.10 

3.60 

2.80 

OB3NV 

10.0 

5,70 

3.36 

2.70 

00 

6.10 

3.60 

2.80 

RB-CO 

1.0 

3.80 

2.10 

2.30 

OD 

6.10 

3.60 

2.70 

RK-CN 

6.0 

5.64 

3.47 

2.70 

00 

6.15 

3.64 

2.80 

SZ-NV 

2.0 

6.10 

3.60 

2.70 

00 

6.10 

3.60 

2.70 

VF-NV 

0.18 

1.30 

0.659 

1.75 

0.55 

2.00 

1.07 

1.196 

00 

5.70 

3.36 

2.70 

YF2NV 

0.24 

1.30 

0.659 

1.75 

0.58 

2.00 

1.07 

1.196 

00 

5.70 

3.36 

2.70 

YF3NV 

0.29 

1.30 

0.659 

1.75 

0.61 

2.00 

1.07 

1.196 

a 

5.70 

3.36 

2.70 

YF4NV 

0.29 

1.30 

0.659 

1.75 

0.70 

2.00 

1.07 

1.196 

Ob 

5.70 

3.36 

2.70 
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1977),  which  was  assumed  to  have  the  same  elastic  properties  as  the  granite 
stock  at  0B2NV  (which  was  modeled  by  a  simple  homogeneous  half space).  This 
similarity  makes  the  pulse  sizes  of  the  computed  synthetic  records  comparable 
with  correction  for  halfspace  properties  below  the  layered  stack. 

Finite  difference  methods  were  also  used  to  estimate  crustal  amplifica¬ 
tion  at  the  four  stations  on  Yucca  Flats  (the  method  is  described  in  the 
paper  by  Kelly  et  al,  1976).  We  attempted  to  model  Yucca  Flats  with  a 
structure  derived  from  Hays  and  Murphy  (1971),  utilizing  Ramspott's  and 
Howard's  (1975)  and  Fernald  et  al's  (1968)  velocity  and  structural  data. 

The  results  of  both  methods  are  shown  in  Table  IV.  Both  the  Haskell  matrix 
and  the  finite  difference  methods  suggest  considerable  amplification  at 
Yucca  Flats,  and  the  two  methods  yield  about  the  same  result. 

Plots  of  m^  differentials  relative  to  0B2NV  corrected  for  crustal 
effects.  Figure  8,  show  major  changes  compared  to  the  uncorrected  data  of 
Figure  6.  The  differences  between  0B2NV  and  the  Yucca  Flats  (YF)  and 
Pahute  Mesa  (NT)  stations  are  drastically  reduced  with  the  exception  of  YF3 
and  YF4  to  a  level  below  signif icance.  The  m^  values  for  FANV,  RBCO  and 
GBNV  are  reduced  below  the  level  of  0B2NV.  The  thin  (h  <  50  m)  weathered 
layer  at  HNME  does  not  justify  any  crustal  correction  relative  to  IFME 
since,  in  our  experience,  such  thin  layers  do  not  amplify  teleseismic  waves. 
However,  any  such  correction  could  move  HNME  to  slightly  lower  m^  values. 

The  relative  positions  of  RKON  and  0B2NV,  the  two  key  stations  on  granite, 
are  unchanged,  and  the  differential  in  m^,  0.17  m.u.,  is  statistically 
significant  -  at  higher  than  the  99%  confidence  level.. 

The  crust-corrected  trace  amplitudes  expressed  in  base  10  log  units 
behave  similarly  as  shown  in  Figure  9.  The  crustal  corrections  for  trace 
amplitudes  are  the  same  as  for  m^  since  there  is  no  apparent  change  in 
period  in  the  synthetic  calculations.  The  0B2NV-RK0N  differential  of  0.267 
is  highly  significant  statistically  and  is  in  excellent  agreement  with  the 
results  of  Butler  (1979). 
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TABLE  IV 


Results  of  finite  difference  and  Haskell  algorithm  calculations  of  crustal 

response  at  Yucca  Flats 

Stations  Magnitude  Differential  Relative  to  0B2NV 


Observed 

Layered  (Haskell) 

FINITE  DIFFERENCE 
1  Hz 

YF 

0.306 

0.34 

0.26 

YF2 

0.370 

0.34 

0.32 

YF3 

0.434 

0.32 

0.27 

YF4 

0.413 

0.32 

0.30 

-31- 


Possible  Bias  Effects  in  m^  Due  To  Variable  Thresholds,  Noise  Level  and 
Source  Region 

Since  relative  and  absolute  amplitude  measurements  can  be  biased  by 
threshold  and  source  region  effects  (von  Seggern  and  Blandford,  1972;  Ring- 
dal,  1976),  we  have  tested  our  data  to  determine  the  extent  of  any  possible 
bias.  The  criteria  we  have  used  for  amplitude  measurement  require  that  S/N 
be  greater  than  3.  Thus,  event  pairs  were  eliminated  when  one  of  the  stations, 
e.g.  station  A,  has  S/N  <  3.  Because  this  condition  could  result  from  low 
signal  at  station  A  relative  to  station  B,  the  average  amplitude  at  A  could 
be  overestimated  relative  to  B  at  low  magnitudes.  Only  if  the  great  majority 
of  all  readings  are  considerably  above  the  3:1  S/N  level  can  such  bias  in 
the  average  be  ruled  out.  To  test  for  this  bias,  we  plotted  the  differen¬ 
tials  of  for  selected  stations  pairs  against  the  average 
two  stations.  Pronounced  trends  in  plots  such  as  these  would  indicate  biases 
in  the  procedures  used  for  determining  An^.  A  slight  change  of  An^  with 
increasing  magnitude  might  also  suggest  a  shift  of  corner  frequency  of 
seismic  sources  of  lower  frequencies.  This  shift  would  be  more  visible  at  a 
high  Q  than  at  a  low  Q  station.  Absence  of  a  clear  trend  indicates  that 
bias  in  Am^  from  variable  noise  levels  is  not  significant. 

Figures  10  through  16  show  such  plots  for  a  selected  set  of  key  station 
pairs.  None  of  these  plots,  including  the  critical  pair  RK0N-0B2NV,  suggest 
a  clear  trend.  Therefore,  bias  effects  from  our  procedure  are  probably 
negligible.  Note  that  the  noise  level  is  approximately  0.3  m.  u.  higher  at  RKON 
than  at  0B2NV.  However,  the  raw  amplitudes  on  the  film  (magnitude  at  1  Hz, 
not  amplitude  corrected  at  T  or  at  A/T)  also  average  0.3  m.  u.  at  RKON.  Thus, 
the  average  expected  S/N  is  the  same,  and  no  bias  would  be  expected. 

To  test  for  the  effect  of  source  region  biases  on  the  measured  relative 
m^,  we  divided  the  RK0N-0B2NV  interstation  differentials  into  four  groups 
marked  with  different  symbols  as  follows: 

Symbol 

A  (Am  -  u)  <  -  a  high  amplitude  at  0B2NV 

b 

B  -  a  <  (Am^  -  p)  <  0 

C  0  <  (Am^  -  p)  <  +  a 

D  (Am,  -  p)  >  +  a  low  amplitude  at  0B2NV 


ror  the  same 
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RK-ON 


Interstation  m,  differentials  plotted  against  the  two 
station  average  m,  for  station  pair  RKON  -  0B2NV.  Th 
absence  of  any  clear  increasing  or  decreasing  trend 
indicates  that  magnitude  threshold  biases  are  not 
present  (see  text). 


GB-NM  0B2NV 


♦  ♦  +  ♦  + 
++  %  +  ♦♦ 


+  %♦  *£+*>+  * 
+  •*+  *  +  * 

+  *  +  * 

♦  ♦  <*++  + 

♦  ♦  V  ♦♦ 


+  + 


+ 

♦ 


+  + 


1.0 


4.0 


5.0 


OB2 


mb 


2 


6.0 


Figure  13  Inter station  m,  differentials  plotted  against  the  two- 

station  average  for  the  station  pair  GBNM  -  0B2NV. 

The  absence  of  any  clear  Increasing  or  decreasing  trend 
indicates  that  magnitude  threshold  biases  are  not  present 
(see  text). 


2 


Figure  14  Interstation  differentials  plotted  against  the  two- 
station  average  for  the  station  pair  NT2NV  -  0B2NV. 

The  absence  of  any  clear  Increasing  or  decreasing  trend 
indicates  that  magnitude  threshold  biases  are  not  present 
(see  text). 


-39- 


NT-NV  0B2NV 


Figure  16  Interstation  a.  differentials  plotted  against  ttewo- 
station  average  a.  for  the  station  pair  NTNV  -  MS*. 

The  absence  of  at#  clear  increasing  or  decreasing  trend 
indicates  that  magnitude  threshold  biases  are  not  present 

(see  text). 


and  plotted  these  symbols  against  azimuth  and  distance  to  0B2NV.  Any  re¬ 
gional  bias  such  as  a  consistent  orientation  of  fault  planes  would  cause 
similar  symbols  to  cluster  in  some  regions.  No  such  clustering  is  discern¬ 
ible  on  the  diagram  shown  in  Figure  17,  and  therefore  there  are  no  indications 
of  source  related  biases. 

Elimination  of  source  region  bias  does  not  mean  that  near-station 
azimuthally  dependent  focusing  on  small  scale  crustal  inhomogeneities  is  not 
present  in  the  data.  Such  biases  cannot  at  present  be  removed  from  these 
data  sets  since  most  of  the  events  occur  in  belts  of  seismicity,  and  thus 
the  azimuth-distance  distribution  of  events  is  uneven. 


Anelastic  Attenuation  as  the  Cause  of  the  RK0N-0B2NV  Magnitude  Differential 

Having  demonstrated  the  existence  of  magnitude  differentials  at  the 
various  SDCS  stations,  we  turn  now  to  their  cause.  In  our  opinion,  anelastic 
attenuation  is  responsible  for  most  if  not  all  of  the  observed  differentials. 
However,  since  station  amplitude  or  magnitude  anomalies  are  subjected  to  biases 
by  such  factors  as  crustal  inhomogeneities,  it  is  important  to  find  diagnostics 
that  indicate  that  the  magnitude  differentials  are  indeed  due  to  anelastic 
attenuation. 

Concentrating  on  the  critical  RKGN-0B2NV  pair  (often  used  in  key  argu¬ 
ments  in  this  report),  strong  support  for  anelastic  attenuation  can  be  found 
by  the  simple  technique  of  subdividing  the  event  population  at  RKON  into  high 
frequency  (broad  band)  signals  and  low  frequency  (narrow  band)  signals  based 
on  the  dominant  period.  The  two  populations  are  identified  as  T  <  T  and 
T  >  T,  where  T  is  the  average  period  for  the  station  for  all  events.  These 
populations  should  behave  as  predicted  in  the  discussion  in  section  B  of 
Part  III  of  this  report.  This  appears  to  be  the  case.  As  shown  in  Figure  18, 
the  histogram  of  base  10  logarithms  of  trace  amplitude  differentials  corrected 
for  distance  shows  that  the  differential  in  trace  amplitude  is  greater  for  the 
population  T  <_  T,  which  agrees  with  the  hypothesis  that  high  frequency  events 
undergo  greater  attenuation.  Histograms  of  Am^  (magnitudes  computed  without 
dividing  by  the  period)  show  reduced  differentials  for  both  populations,  and 
at  the  same  time  the  difference  in  the  means  of  the  two  populations  also 
decreases.  This  is  due  to  the  over-correction  associated  with  the  assumption 
of  a  monochromatic  signal  (see  section  B).  Finally,  the  last  line  of 


Figure  17  Plot  of  magnitudes  of  interstation  a  differentials  for 

the  RKDN-0B2NV  pair  against  distance  and  azimuth  of  0B2NV . 
The  absence  of  any  clear  clustering  of  symbols  in  any 
region  indicates  that  there  are  no  source  region  biasing 
effects  due  to  dominant  source  orientations  (see  text). 


Am.'  A  (*09,0  Atr) 


Figure  18  Subdivision  of  the  event  population  into  wide  band  and 
narrow  band  signal  populations  shows  that  a)  the  RKON- 
0B2NV  trace  amplitude  differential  is  less  for  narrow 
band  longer  period  signals  (top  row),  b)  the  differential 
is  reduced  in  m',  trace  amplitude  divided  by  the  instru¬ 
ment  magnification  at  the  dominant  period  measured.  This 
reduction  is  greater  for  wide  band  signals  (middle  row) 
and  c)  division  by  period  T  causes  the  increase  of  m^ 
differential  relative  to  that  of  m' ,  and  the  increase  is 
greater  for  wide  band  signals.  This  behavior  is  diagnos¬ 
tic  of  attenuation  as  a  cause  for  the  RK0N-0B2NV  magni¬ 
tude  differential  (see  text  and  section  B  of  Fart  III). 
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Figure  18  shows  the  computed  by  dividing  by  the  period.  This  tends  to 
restore  the  regional  differential  since  the  period  is  shorter  at  RKON  than 
at  0B2NV.  The  period  differential  is  also  greater  for  the  population  T  £  T 
as  predicted  in  section  B  and  Judging  by  the  increase  of  Am^  versus  im\ 

This  behavior  is  quite  consistent  with  the  interpretation  of  the 
0B2NV-RK0N  magnitude  differential  as  an  attenuation  effect,  and  it  effec¬ 
tively  rules  out  many  other  interpretations  such  as  accidental  local  focusing. 
One  feature  of  this  reasoning  is  that  it  is  based  entirely  on  time  domain 
measurements  and  should  suit  those  who  dislike  spectral  arguments. 
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SPECTRAL  ANALYSIS 


Relative  t*  values  (At*)  for  the  SDCS  stations  have  been  computed  by 
the  spectral  ratio  method.  Because  computing  spectral  ratios  for  all 
possible  pairs  of  stations  is  redundant,  we  computed  ratios  for  station 
pairs  directly  connected  by  lines  in  Figure  19.  Differentials  in  t*  and 
their  standard  deviations  for  other  station  pairs  can  be  derived 

easily  from  these  results.  This  approach  takes  advantage  of  the  fact  that 
for  closely  located  station  pairs,  has  smaller  variance  than  for  distant 

pairs  of  stations.  (This  will  be  discussed  in  more  detail  later.)  Thus, 
for  example,  determining  At*  for  the  RK0N-0B2NV  and  0B2NV-YFNV  pairs  is 
more  reliable  than  determining  At*  directly  from  the  few  events  common  to 
RKQN  and  YFNV  alone.  This  approach  is  also  the  most  practical  because 
during  the  project's  two  phases  many  events  for  the  RK0N-0B2NV  pair  were 
available  to  reduce  the  variance  of  their  mean  At*,  while  only  a  few  events 
were  sufficient  to  define  At*  for  the  0B2NV-YFNV  pair  to  the  same  accuracy. 

The  histograms  of  the  measured  At*  are  given  in  Appendix  D.  In  addi¬ 
tion  to  those  involving  the  new  stations  of  Phase  III,  updated  versions  of 
histograms  for  the  0B2NV-RKQN  and  HNME-RKQN  pairs  are  presented.  In  all 
cases,  an  accuracy  (2a^t*)  of  'UJ.OS  sec  was  the  goal. 

The  easiest  way  to  discuss  the  At*  values  is  to  compare  them  to  a  common 
standard  station.  Figure  20  summarizes  the  results  using  0B2NV  as  the 
standard.  As  mentioned  above,  some  of  these  At*  mean  values  and  their 
standard  deviations  were  derived  indirectly.  This  figure  also  includes  some 
relative  t*  data  at  SZNV  (SHOAL  site,  on  granite)  and  SEMN  with  the  assump¬ 
tion  that  t*  at  SEMN  is  the  same  as  at  RKON.  This  can  be  justified  easily 
by  S  wave  data  shown  later  in  this  report.  Figure  20  shows  that  all  WUS 
stations  have  essentially  the  same  t*  as  0B2NV  except  stations  NTNV  and 
NT2NV,  which  have  slightly  lower  t*  than  0B2NV.  On  the  other  hand,  the 
RK0N-0B2NV  differential  in  t*  is  about  0.2  sec  and  highly  significant  statis¬ 
tically.  The  HNME-0B2NV  differential  is  less  (0.1  sec),  but  it  is  also 
significant  at  the  95%  confidence  level.  This  lower  value  indicates  some 
attenuation  under  the  northeastern  U.S.,  a  possibility  suggested  by  Solomon 
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to  that  of  0B2NV.  RKON  on  a  shield  has  the  lowest  t* 
while  the  WUS  stations  have  the  largest.  95%  confi¬ 
dence  limits  are  indicated  by  bars. 


and  Toksoz  (1970)  and  Der  et  al  (1975).  The  t*  at  IFME  is  comparable  to  that 
at  HNME. 

There  was  some  difficulty  with  the  Yucca  Flat  stations  in  reliably  deter¬ 
mining  t*.  This  difficulty  resulted  from  high  unstationary  noise  levels  due 
to  drilling.  As  a  result,  the  standard  procedure  for  estimating  noise  in  the 
signal  window  (taking  spectra  of  time  windows  prior  to  the  arrival  of  the  P 
wave)  failed  repeatedly.  Many  of  the  spectral  ratios  supposedly  show  points 
with  high  S/N  (around  4  Hz)  that  disagree  with  the  trend  at  lower  frequencies. 
This  can  be  attributed  to  nonstationary  noise  rather  than  to  real  signal 
energy  at  4  Hz,  which  should  be  far  beyond  the  comer  frequency  expected  for 
most  events.  Moreover,  crustal  response  calculations  rule  out  major  enhance¬ 
ments  of  high,  frequency  energy.  Thus,  the  reader  should  assign  less  weight 
to  the  t*  determination  at  YF  stations. 

Returning  to  Figure  20  the  OB2NV-OB3NV  differential,  although  small,  is 

significant  statistically,  t*  for  stations  NTNV,  NT2NV  and  GQNV  appears  to 

be  smaller  than  that  at  0B2NV  and  Yucca  Flats  although  significantly  higher 

than  that  of  RKON .  A  high-Q  high-velocity  plug  may  be  present  under  Pahute 

Mesa  as  suggested  by  Spence  (1974),  and  this  could  also  affect  measurements 

at  nearby  GQNV.  The  high  frequency  content  of  L  is  also  visibly  enhanced  at 

8 

Pahute  Mesa  (Barker  et  al,  1980).  Since  propagates  entirely  in  the  crust, 
this  may  indicate  that  the  apparently  low  t*  at  NTNV  and  NT2NV  could  be  no 
more  than  a  local  crustal  resonance  effect.  We  have  no  preferred  interpreta¬ 
tion  of  this. 

As  shown  in  Figure  21,  the  dominant  periods  of  P  waves  behave  similarly 

to  the  At*  values  of  Figure  20.  AT  is  smaller  at  NTNV,  NT2NV,  and  GQNV  than 
P 

at  0B2NV  and  the  Yucca  Flats  stations.  Moreover,  all  of  the  WUS  stations 
have  AT  values  significantly  higher  than  that  of  RKON.  Periods  at  HNME  are 
not  comparable  to  those  of  other  SDCS  stations  due  to  the  difference  in 
instrumentation  at  HNME.  Histograms  of  AT  for  the  rest  of  the  SDCS  Network 
are  compiled  in  Appendix  E. 

As  mentioned  earlier,  the  standard  deviation  of  relative  t*  values  also 

depends  upon  mutual  distance  between  stations  in  a  manner  similar  to  that 

shown  by  the  magnitude  residuals.  A  plot  of  o  ^  versus  A0  is  shown  in 

Figure  22.  Note  that  the  numerical  values  of  a  are  smaller  and  increase 

At* 
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DISTANCE  DEGREES 


more  slowly  with  distance  than  oA  «  Because  absolute  numerical  values  of 
Am^  and  At*  are  both  of  the  order  of  a  few  tenths  in  magnitude  units  and 
seconds,  respectively,  fewer  measurements  are  needed  to  establish  At* 
between  stations  than  to  determine  magnitude  residuals.  This  result  is  in 
agreement  with  experience  at  seismic  arrays  and  shows  that  At*  is  numeri¬ 
cally  a  more  stable  quantity  than  Am^  in  terms  of  multipathing  and  other 
disturbances  (Der  et  al,  1977). 

Narrow  Band  Determination  of  t*  Differential  for  RKON  and  0B2NV 

In  order  to  see  whether  the  At*  value  we  computed  depends  excessively 
on  the  low  level  high  frequency  end  of  the  spectrum  (3  to  4  Hz) ,  we  recom¬ 
puted  the  relative  t*  between  RKON  and  0B2NV  using  spectral  slope  fits  in  the 
narrower  0.5  to  2  Hz  band.  A  strong  frequency  dependent  variation  in  the 
relative  t*  differential  between  0B2NV  and  RKON  would  result  in  a  drastically 
different  t*  differential  for  this  narrower  band.  Figure  23  shows  that  this 
is  not  the  case — the  relative  t*  differential  0.24  ±  .06  sec  is  not  signifi¬ 
cantly  different  from  the  wider  band  value  of  0.20  ±  .03  sec.  This  demon¬ 
strates  that  our  results  do  not  depend  on  low  level  high  frequency  energy, 
the  existence  of  which  is  questioned  by  some  researchers.  This  statement  is 
also  supported  by  the  fact  that  time  domain  measurements  that  depend  primarily 
on  the  1  to  2  Hz  band  are  also  indicative  of  t*,  variation  between  RKON  and 
0B2NV.  The  above  result  also  rules  out  a  rapid  change  with  frequency  in  the 
Interstation  t*  differential  above  1  Hz. 

Testing  for  some  Biases  in  Relative  t*p  Measurements 

To  see  whether  any  source  region  bias  is  apparent  in  the  measured  rela¬ 
tive  t*,  we  separated  the  RK0N-0B2NV  At*  values  into  four  groups  designated 

by  different  symbols  as  follows: 

Symbol 

A  (At*  -  y)  <  -  o 

B  -a  <  (At*  -  y)  <  0 

C  0  <  (At*  -  y)  <  +  a 
D  (At*  -  y)  >  +  o 
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Figure  23  Histograms  of  0B2NV-RK0N  t*  differentials  computed  in 

the  0.5  to  2.0  Hz  and  the  0.5  to  4.0  Hz  bands.  The  fact 
that  the  narrow- band  average  differential  is  about  the 
same  as  that  computed  in  the  0.5  to  4.0  Hz  range  (0.244 
versus  0.200)  rules  out  rapid  change  of  relative  t*  with 
frequency  in  this  range,  and  shows  that  our  t*  results 
do  not  depend  critically  on  low  level  high  frequency 
energy.  Note  that  the  scatter  in  these  t*  histograms  is 
small  compared  to  that  in  Figures  3  and  4  for  m^  and  trace 
amplitudes,  demonstrating  the  greater  consistency  of 
spectral  measurements. 
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and  plotted  these  symbols  against  azimuth  and  distance  to  0B2NV  In  Figure  24. 
This  subdivision  assigns  the  letter  D  to  the  largest  t*  difference  and  A  to 
the  least.  The  resulting  pattern  shows  no  convincing  evidence  of  regional 
biases.  Such  biases  could  in  theory  be  due  to  directionality  of  finite 
sources  that  would  result  in  consistent  regional  patterns  of  At*  if  the 
fault  planes  were  oriented  in  any  consistent  manner  relative  to  the  station 
pair  studied.  The  lack  of  any  such  bias  effect  is  not  surprising  since, 
although  the  geographical  distribution  of  sources  used  is  clearly  non-random, 
the  orientation  of  the  trends  of  seismic  belts,  island  arcs  and  the  like 
relative  to  the  0B2NV-RK0N  station  pairs  varies  enough  to  eliminate  such  bias 
even  if  the  source  mechanism  of  all  events  were  the  same.  But  the  fact  is 
that  sources  in  the  seismic  belt  include  a  variety  of  mechanisms,  and  this 
further  reduces  chances  for  finding  overall,  consistent  source  related 
biases  in  the  determination  of  At*  from  averages  of  many  events. 

To  demonstrate  further  that  At*  behaves  consistently,  consider  the 
recordings  of  a  Novaya  Zemlya  explosion  in  Figure  25.  On  the  left  is  the 
recording  at  0B2NV  and  on  the  right  at  RKQN.  This  source  is  supposedly 
symmetrical  azimuthally,  but  the  trace  amplitudes  differ  by  a  factor  of  30, 
a  dramatic  illustration  of  the  instability  of  amplitude  measurements.  The 
filtered  traces  shown  below  demonstrate  that  the  P  wave  is  richer  in  high 
frequencies  at  RKON  than  at  0B2NV,  The  S/N  ratio  is  high  for  this  event. 

This  figure  also  demonstrates  that  the  concept  of  "transparency"  used  in 
time  domain  work  is  not  a  physically  meaningful  one.  0B2NV  was  regarded  as 
a  "transparent"  station  in  much  of  the  time  domain  work  of  Butler  et  al 
(1979)  and  it  is  clearly  not  in  this  case.  The  waveform  of  P  is  quite 
complex.  Experience  at  seismic  arrays  indicates  that  "transparency"  is  a 
function  of  near  receiver  focusing  and  thus  varies  with  the  azimuth  of  the 
events  observed. 

Studies  of  S  Waves  at  SDCS  Stations 

Although  the  prime  emphasis  in  this  project  was  on  P  wave  amplitude 
residuals  and  spectra,  a  limited  search  was  made  for  short-period  S  waves 
recorded  in  the  SDCS  data  base.  The  S  wave  data  thus  found  is  limited,  and 
the  S/N  ratio  of  these  signals  is  generally  not  good.  In  Figures  26  through 
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31  we  show  a  few  examples  of  short-period  S  waves  and  the  spectral  ratios 
computed  from  them,  p-  to  the  poor  S/N  ratios  of  these  signals  most  of  the 
spectral  ratios  arc  __ jd  over  a  frequency  range  of  only  0.75  Hz.  Neverthe¬ 
less,  the  available  data  indicate  that  RKON  is  richer  in  high  frequency  con¬ 
tent,  and  GBNM  appears  to  have  t*  values  comparable  to  0B2NV.  Due  to  the 

s 

sparse  data  and  the  emphasis  on  m^,  S  wave  amplitude  residuals  were  not 
studied  at  SDCS  stations. 
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ratio  for  the  station  pair 
iponent),  19  June  1977 


S  wave  spectra  and  spectral  ratio  for  the  station  pair 
RKON  -  0B2NV  (transverse  component),  19  June  1977 
11:47:22.3,  Kurile  Islands. 


TRAVEL  TIME  RESIDUALS 


Since  P  wave  travel  time  residuals  constitute  an  additional  parameter 
that  characterizes  the  physical  state  of  the  mantle  under  a  given  station, 
we  have  computed  relative  travel  time  residuals  for  all  SDCS  stations  with 
sufficient  data  base.  Late  arrivals  are  commonly  associated  with  high  upper 
mantle  attenuation  and  the  presence  of  a  low  velocity  layer  in  the  upper 
mantle.  Worldwide  observations  support  the  correlation  of  travel  time  resi¬ 
duals  with  highly  attenuative  properties  and  other  diagnostic  geophysical 
measurements  such  as  S  wave  delays,  heat  flow,  and  electrical  conductivity. 
In  view  of  this  fact,  it  is  important  to  measure  the  relative  travel  time 
residuals  among  the  SDCS  stations  in  order  to  further  evaluate  the  condition 
of  the  upper  mantle  under  each  station.  Travel  times  for  P  waves  were 
routinely  compiled  throughout  the  SDCS  project,  but  it  is  only  recently 
that  we  have  accumulated  sufficient  data  to  compute  travel  time  residuals. 

The  stations  of  prime  interest  are  0B2NV,  RKON  and  IFME,  all  of  which 
are  located  on  granite,  though  great  distances  apart.  Determination  of 
relative  residuals  between  the  members  of  this  group  is  complicated  by  large 
scatter  in  the  results  of  the  computation.  Generally  speaking,  the  greater 
the  distance  between  stations  the  greater  the  scatter.  Thus,  the  problem 
is  less  serious  for  the  group  of  stations  located  throughout  the  WUS 
(0B2NV,  FANV,  GBNM)  since  the  relative  distances  are  less.  Finally,  the 
stations  at  NTS  constitute  a  tightly  spaced  group  for  which  the  measurement 
of  relative  residuals  constitutes  no  difficulty. 

In  all  travel  time  studies  one  also  has  to  deal  with  some  measurements 
which  are  grossly  in  error.  This  did  not  prove  to  be  extremely  difficult. 

In  most  cases,  as  a  large  number  of  readings  was  compiled,  we  found  a  tight 
group  resembling  a  normal  distribution  with  some  very  obvious  outliers.  We 
followed  the  standard  statistical  practice  of  computing  the  standard  devia¬ 
tion  a  for  the  tight  group  and  omitting  all  points  removed  by  mote  than  3cr 
from  this  group.  The  mean  travel  time  residuals  were  then  computed  from  the 
purged  data  set.  Histograms  of  travel  time  differences  for  various  station 
pairs  are  given  in  Appendix  F.  The  differentials  and  their  95%  confidence 
limites  are  tabulated  in  Table  V. 


TABLE  V 

Travel  time  residuals  for  selected  SDCS  stations  relative  to  0B2NV 


STATION 

ATT 

• 

ATTC°rr 

FANV 

0.03 

+ 

0.16 

-0.52 

GBNM 

-0.05 

± 

0.25 

-0.49 

HNME 

-1.28 

± 

0.63 

-1.04 

NT  NV 

0.43 

± 

0.23 

-0.23 

RKON 

-1.89 

+ 

0.29 

-1.68 

YFNV 

0.61 

+ 

0.22 

0.37 
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As  mentioned  above,  the  standard  deviations  of  travel  time  differentials 
have  a  distance  dependence  similar  to  those  of  the  magnitude  and  t*  differen¬ 
tials.  This  is  shown  in  Figure  32.  The  physical  causes  for  this  phenomenon 
are  the  increasing  effect  of  mislocation  of  events  with  distance  and  the  in¬ 
creasing  probability  of  misreading  the  first  arrivals  due  to  the  decreasing 
visual  similarity  of  the  signals. 

Table  V  also  shows  the  residuals  corrected  for  elevation  and  local 
geology  and  adjusted  relative  to  0B2NV  (0B2NV  is  assumed  to  be  zero).  The 
table  shows  that  all  WUS  stations  are  late  relative  to  RKON,  which  agrees  with 
results  obtained  by  other  workers  (Sengupta  and  Julian,  1976).  Stations  in 
the  northeastern  United  States,  HNME  and  IFME,  are  later  than  RKON  but  not  by 
as  much  as  the  WUS  stations.  This  picture  is  in  perfect  agreement  with  other 
studies  and  also  conforms  to  our  evaluation  of  relative  t*  differentials, 
assuming  that  LVZ  is  associated  with  low  Q. 
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SUMMARY 

Numerical  values  of  the  most  relevant  diagnostic  parameters  obtained  In 
the  data  analysis  are  summarized  In  Table  II.  The  data  Indicate  that  after 
crustal  corrections  are  made,  P  wave  amplitudes  tend  to  be  lower  at  the  WUS 
test  sites  than  at  RKON.  The  NEUS  stations,  HNME  and  IFME,  appear  to  be 
intermediate  in  amplitude  between  the  two  groups,  but  large  scatter  in  the 
data  does  not  permit  the  determination  of  a  significant  m^  difference  between 
these  stations  and  either  0B2NV  or  RKON.  Crustal  corrections  leave  large 
differences  between  the  stations  in  or  trace  amplitudes  which  are  unex¬ 
plained  in  terms  of  t*  or  known  structures  and  are  presumably  caused  by  local 
focusing  effects. 

Spectral  measurements  indicate  a  significant  loss  of  high  frequency  con¬ 
tent  at  all  WUS  test  sites  relative  to  RKON.  The  NEUS  stations,  HNME  and 
IFME,  occupy  an  intermediate  position  in  frequency  content.  The  t*  deter¬ 
mined  from  spectral  ratios  shows  a  good  qualitative  agreement  with  the 
dominant  periods  measured  by  analysts. 

Travel  time  residuals  shown  in  Table  V  indicate  that  after  local  crustal 
and  elevation  corrections,  P  waves  arrive  late  at  all  WUS  test  sites  relative 
to  RKON.  The  NEUS  stations,  HNME  and  IFME,  again  appear  to  be  intermediate 
in  their  relationship. 

The  signific  -,nce  of  these  results  will  be  discussed  in  the  next  part  of 
this  report  in  the  broader  context  of  related  work  by  SDAC  and  other 
researchers. 


-68- 


PAST  2 

EVALUATION  OF  THE  RESULTS  OF  THE  SDCS  PROJECT  IN  THE  CONTEXT  OF  RELATED  WORK 

BY  SDAC  AND  OTHER  RESEARCHERS 


The  first  part  of  this  report  presented  the  results  of  data  analysis 
from  the  SDCS  and  a  few  LRSM  stations.  The  areal  coverage  of  this  station 
set  is  limited;  only  one  of  the  stations,  RKQN,  is  located  on  a  shield.  Thus 
it  is  extremely  important  to  show  that  RKQN  is  not  in  any  way  anomalous  lest 
our  conclusion  concerning  the  other  shield  regions  be  seriously  in  error. 

The  station  in  Maine  should  also  be  tested  in  the  same  context.  Since  there 
have  been  other  studies  of  amplitude  and  magnitude  anomalies  of  teleseismic 
P  waves,  our  results  must  be  interpreted  in  the  framework  provided  by  these. 
For  example,  some  recent  broad  regional  studies  of  P  wave  anomalies  were 
interpreted  in  a  manner  that  seemed  to  be  in  conflict  with  our  conclusions. 

We  shall  show  that  no  real  conflict  exists  if  some  errors  and  unintentional 
misrepresentations  are  corrected.  We  also  include  here  some  results  of  our 
own  spectral  and  amplitude  studies  of  P  and  S  waves.  In  addition,  we  shall 
outline  regional  variations  of  under  the  contiguous  United  States.  In 
the  last  part  of  this  section  we  shall  put  some  limits  on  the  absolute  values 
of  t*  and  t*  and  discuss  possible  forms  of  frequency  dependence  of  t*  for 
two  types  of  paths. 


Discussion  of  the  Results  of  the  SDCS  Experiment  in  the  Context  of 
and  Spectral  Studies  of  Short-Period  P  and  S  Waves 


ilitude 


Although  we  pointed  out  In  the  previous  part  of  this  report  that  the 
SDCS  data  can  be  interpreted  in  terms  of  decreased  Q  in  the  upper  mantle, 

p 

there  is  still  a  need  to  integrate  these  results  into  a  framework  of  other 
regional  studies  to  confirm  our  interpretation.  Any  lateral  decrease  of  Q 

P 

in  the  mantle  would  cause  the  following  suite  of  phenomena  all  of  which  have 
to  be  present: 

a)  Decrease  of  P  wave  amplitudes  resulting  in  a  magnitude  anomaly 
for  P  waves. 

b)  Decrease  of  high  frequency  content  in  P  waves. 

c)  A  regional  anomaly  for  S  wave  amplitudes.  The  effect  should 
be  greater  than  that  for  P  waves. 

d)  High  frequency  content  of  S  waves  should  decrease  at  a  greater 

degree  than  that  for  P  waves  since  t*  -  4t*. 

s  p 
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We  shall  discuss  manifestations  of  these  phenomena  and  show  that  all  of  these 
are  present  in  much  of  the  western  United  States. 

P-Wave  Amplitude  Anomalies 

It  would  be  expected  that  in  the  regions  underlain  by  a  low  Q  mantle  the 
P-wave  amplitudes  from  observed  teleseisms  would  be  reduced.  This  has  been 
confirmed  in  several  studies  of  magnitude  variations  across  the  U.S.  (Evernden 
and  Clark,  1970;  Cleary,  1967;  North,  1977;  Booth  et  al,  1974).  The  various 
studies  all  show  a  negative  bias  in  amplitudes  in  most  of  the  WUS  relative  to 
the  EUS  as  a  whole.  The  methods  of  data  selection  and  amplitude  measurement 
methods  vary,  and  the  reported  bias  values  for  common  stations  in  many  cases 
also  differ  considerably. 

Figure  33  shows  the  regional  pattern  obtained  by  Booth  et  al  (1974). 

This  map  shows  strong  negative  magnitude  residuals  in  the  southwestern  United 
States  and  positive  residuals  in  the  north  central  section  of  the  country.  A 
study  of  these  residuals  (Der  et  al,  1979)  showed  that  crustal  amplification 
alone  cannot  explain  this  pattern,  and  that  even  after  correcting  for  crustal 
effects,  a  0.33  m. u.  EUS-WUS  regional  bias  remains.  Figure  34  from  this 

c 

study  shows  that  magnitude  residuals  plotted  against  logarithms  Am^  of 

multiplicative  crustal  amplification  factors  A£  tend  to  cluster  around  two 

regression  lines  for  the  WUS  and  EUS  populations  respectively.  Therefore, 

a  multiple  regression  using  both  the  crustal  amplification  factor  Ac  and 

t*  is  necessary  to  reduce  the  variance.  The  regression  on  the  set  of  sta- 
P 

tions  gave  for  an  assumed  EUS-WUS  t*  difference 

Ara^  "  (1.35  ±  0.32)  At* 

at  the  95%  confidence  level. 

This  coefficient  is  of  the  same  order  as 

nf/lnlO  *  1.36  (for  f  =  1  Hz) 

which  would  be  the  multiplicative  coefficient  between  t*  and  the  10  base 
logarithms  of  amplitudes  for  an  attenuated  1  Hz  wave.  This  equation  shows 
that  average  regional  variations  of  amplitudes  in  short-period  P  and  t*  are 
closely  tied  together  and  cannot  be  specified  independently. 

Consider  now  some  results  of  a  study  by  Butler  and  Ruff  (1980)  in  the 
same  context.  Figures  35  to  39  show  logarithmic  plots  of  P  wave  amplitudes 
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Magnitude  residuals  for  LRSM  stations  (after  Booth  et 
al,  1974). 


IN  KURILE  IS.,  JAPAN,  BONIN  IS. 


TJ  C 

w  a  c 

<U  -H  *H 

W 

T3 

• 

a  q 

X2  4-i 

C 

w 

CO  B 

4J  tI  4J 

C0 

4-i 

B 

U  40 

c 

JZ;  O 

M  CO  00 

W 

cu 

CO  CJ 

O  pL.  *H 

M 

> 

O.  -U  03  CO 
1  (U  CiJ 
CO  MW 
U  QJ 

a)  oJ  4-j  co 
U  H  ^  D 
nJ  -m  re 

4-1  PQ 

(U  L  H  • 

>  Q)  eg 

ctj  4J  B  *j 

>  o  x 
co  a  <u 

pH  (tf  4J 

u-i  w  a)  a) 
o  c  40  a) 
o  4J  w 

W  *H 
B  4J  4h 

JZ  aj  O  w 

AJ  4J  c 

•H  W  W  O 
M  ai  *H 
aJ  co  CX  -u 

W)  o  eg  to 
o  a  JS  4-» 
1-4  co  w  w 


<D  <u 
£  0) 
rH  C 
*-H  *H  M 

w  m  a> 

C  3  4J 

O  w 

♦h  a; 

a  > 

M  O  J2 
O  M  4J 
(XM  M 

o  o 
u  <u  c 

CX  u 
W  CO  M 
•W  <u 

■°  gfi 

S.8  ° 

>  w  -o 
*h  e 
oo  co  m 

4J 

<U  CO  w 
L  TJ  3 
CO  C 
0)  CO 

BlrCH 

C  H  W 
O  W 

4J  •  C 
CO  W  tH 
•U  -U  C 

■5a 

iH  D. 
cd  - 
u  dl  C 
and 
d  hi  a 

O  40  CO 
O  4J  *-) 


AZIMUTH 


a  H03X&SU*  Hinos 


<14  cl  3  C 

CO  »M  cd 

a  c  <u  h-j  u 

cd  Q  Xt  *H  MJ  -H 

a  u  o  s:  u 

z:  g  ni  d)  <u 

cfl  o  m  cm  *h  a 

cfl  o  o  s>  < 

i  ^  ^ 

.  T3  O  tH  4-1 

h  a)  to  cd  3 

O  ^  j3  CO 

14-1  O  O  -O  O  Cfl 

00  4J  c  1H 
id  ov  hi  al  4i  6 
0)  ^  MO 

•C 3  00  CO  O  M 

3  -  3  P->mi 

Mi  .-I  0)  O  O 
•H  cd  M  *H  M  0) 

r-i  n)  Ml  p-  M 

p.  4-1  id  in  d 
a  a)  c  m  t4 

a  M  01  T3  C 

M  0)  _  3 

01  01  4-1  Cfl  C  O 

O  rl  U  P  O  £ 

id  4i  d  »  >  * 

M  3  &.  !Z  v4 

4-1  PQ  00  cd 


oi  m  r-t  •  a)  cd 
>  a)  d  K  3 

d  u  0  4i  id  , 

;*  m-i  o  x  o) 

Cd  G  oi  2  £ 

p,  w  u  4i  C  H 

o 

<44  »  a)  o)  i-4 

o  c  j3  o  Mi  •  • 

o  44  01  Cd  01  CO 

8-_|  44  41  41  d 

Ml  144  0)  O  F* 

d  °  ®  ^  3 

41  41  3  •— 4  CM  3 

•i-(  01  01  o  cd  J 

M  0)  1-4  4J  0)  J3 

cd  cfl  a  4j  <»  oi  mi 

oo  o  cd  <d  cd  2  m 

o  P  Id  41  O  43  5 

j  to  n  #  o  4J  d 


J 


as  E-W  cross  sections  of  sorts  across  the  United  States.  These  represent 
amplitude  patterns  for  Russian  test  sites,  earthquakes  from  the  Kuriles,  all 
earthquakes  for  the  NW  azimuth.  South  American  earthquakes,  and  the  sum  of 
all  of  the  data.  For  convenience  of  presentation  they  spaced  the  stations 
equally  on  the  plots,  a  practice  that  gives  an  enhanced  visual  weight  to 
some  groups  of  stations  that  is  out  of  proportion  to  the  small  areas  they 
occupy  (ie,  the  BLA-SCP-GEO-OGD-WES-HNME  group  and  the  BKS-COR-LON  group). 

We  shall  refer  to  this  fact  later.  Inspecting  these  figures  in  succession 
one  can  see  a  common  pattern  (sketched  in  freehand  with  a  solid  line)  that 
is  clearly  associated  with  the  crust-mantle  structure  under  the  United  States. 
The  wide  variations  from  azimuth  to  azimuth  and  source  region  to  source  region 
are  not  surprising  in  view  of  the  focusing  phenomena  described  in  Section  D 
of  the  last  part  of  this  report.  There  is  an  indication  of  another  detail, 
a  low  amplitude  region  along  the  Rocky  Mountain  front  (dashed  line). 

Another  curious  feature  of  Figures  35  to  39  is  that  on  most  of  these 
plots  the  SDCS  stations  0B2NV,  RKON  and  HNME  are  high  in  amplitude  relative 
to  their  surroundings.  This  arouses  the  suspicion  that  in  Butler's  work, 
when  data  from  the  two  networks,  SDCS  and  WWSSN,  were  assembled  together, 
the  mean  m^  levels  of  the  two  networks  were  not  adjusted  properly.  To  in¬ 
vestigate  this  we  attempted  to  tie  0B2NV  in  with  this  pattern  by  computing 
the  magnitude  differential  at  the  common  station  0B2NV  and  station  ANMO 
(ALQ).  Readings  of  trace  amplitude  differentials  at  corresponding  portions 
of  the  P  wavetrain  for  72  events  resulted  in  the  histogram  of  log.^  (Amplitude) 
shown  in  Figure  40.  The  mean  differential  -0.15  ±  .058  (95%  confidence) 
indicates  that  the  relative  positions  of  WWSSN  and  SDCS  stations  in  the  study 
of  Butler  et  al  (1979)  are  not  correct  since  Butler's  differential  between 
0B2NV  and  ALQ  is  much  larger;  it  is  outside  the  confidence  limits  indicated. 
This  confirms  our  suspicion.  Within  each  network,  the  relative  differences 
in  m^  levels  are,  however,  correct  in  Butler's  study  and  for  the  SDCS  network 
his  results  agree  well  with  ours. 

Figure  35  to  39  were  interpreted  by  Butler  et  al  (1979)  as  showing  that 
no  average  EUS-WUS  regional  difference  in  amplitudes  exists,  and  the  visual 
impression  from  these  figures  appears  to  confirm  this  if  one  averages  the 
levels  in  the  WUS  and  EUS  in  their  graphs.  However,  if  one  considers  the 
relative  areas  covered  by  the  stations  in  the  northeastern  United  States, 
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Figure  40  Histogram  of  trace  amplitude  differentials  of  P  waves 
at  0B2NV  and  ALQ  (ANMO) .  This  differential  is 
significantly  different  from  the  value  given  by 
Butler  et  al,  (1980). 


about  1/20  of  the  total  area,  the  visual  weight  given  to  the  stations  in  the 
NEUS  in  Butler's  figures  (where  they  are  evenly  spaced)  is  excessive.  There¬ 
fore  these  figures  are  highly  deceptive.  If  one  takes  averages  weighted 
with  respect  to  the  areas  covered,  the  conclusion  of  many  researchers  that 
the  EUS  as  a  whole  has  larger  or  amplitude  levels  than  the  WUS  still  stands. 

P  Wave  Spectral  Anomalies 

Due  to  the  sensitivity  of  high  frequency  body  waves  to  variations  of 
t*,  spectral  scopes  are  strong  diagnostics  of  regional  variations  of 
attenuation.  While  body  wave  magnitude  measurements  involve  mostly  1  Hz 
energy,  spectra  of  P  waves  often  show  significant  energy  up  to  4  Hz  that 
is  considerably  more  sensitive  to  variations  of  Q  than  1  Hz  energy.  This 
fact,  together  with  the  greater  stability  of  spectral  slopes  with  respect 
to  lateral  inhomo gene i ties,  makes  spectral  measurements  extremely  important 
in  evaluating  mantle  Q  under  any  location. 

Spectral  measurements  can  be  used  in  various  ways  to  estimate  t*  in  the 
earth.  Spectral  ratios  of  the  observed  P  spectra  to  the  estimated  source 
spectrum  can  be  used  to  estimate  the  absolute  value  of  t*.  In  most  cases, 
absolute  strengths  of  sources  are  not  well  known  for  either  earthquakes  or 
explosions;  however,  the  slopes  of  source  spectra  are  well  enough  specified 
to  put  reasonable  limits  on  t*.  For  explosions,  near  source  measurements 
are  often  available  to  estimate  the  source  spectrum  from  a  reduced  displace¬ 
ment  potential.  In  the  absence  of  near  source  data  it  appears  that  source 
models  such  as  those  of  von  Seggem  and  Blandford  (1972)  or  Mueller  and 
Murphy  (1971)  describe  explosion  spectra  with  sufficient  accuracy  to  permit 
the  computation  of  t*  to  within  0.1  second.  For  earthquake  observations, 

basic  physical  arguments  require  that  the  spectra  fall  off  beyond  a  specific 

-2  -3 

(although  ill  defined)  corner  frequency  at  rates  of  w  or  m  ,  a  condition 
that  immediately  puts  limits  on  the  possible  range  of  absolute  t*. 

In  addition  to  estimating  absolute  t*,  spectral  ratios  of  body  waves 
observed  at  various  stations  for  common  events  can  be  used  to  determine 
regional  variations  of  anelastic  attenuation.  The  ideal  sources  for  such 
studies  are  nuclear  explosions  because,  except  for  possible  strain  release, 
they  are  non-directional  sources  essentially  radiating  the  same  spectrum  in 
all  directions. 
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In  our  previous  papers  (Der  and  McElfresh,  1976,  1977)  we  compiled  t* 
estimates  for  a  group  of  paths  crossing  the  mantle  under  the  United  States, 
mostly  for  nuclear  explosions.  Paths  not  crossing  the  mantle  under  the  WUS 
were  associated  with  signals  of  significant  energy  in  the  3  to  4  Hz  range. 
Paths  crossing  the  mantle  under  the  WUS  showed  a  great  reduction  of  the 
high  frequency  content.  The  t*  for  shield  paths  was  found  to  be  around  0.2 
sec  while  for  paths  crossing  the  upper  mantle  under  the  WUS  the  t*  was 
between  0.4  and  0.5  sec.  None  of  our  data  showed  t*  equal  or  larger  than 
1  sec,  although  it  is  possible  that  some  paths  from  the  WUS  to  other  tectonic 
regions  might  have  t*  ^  1  sec.  Spectral  ratios  computed  at  the  SDCS  stations 
also  fall  into  the  pattern  outlined  above.  All  WUS  test  sites  show  relatively 
high  t*  compared  to  RKON,  while  the  stations  in  the  northeastern  U.S.  are 
intermediate  with  regard  to  t*,  a  position  that  correlates  with  the  regional 
amplitude  patterns. 

S  Wave  Amplitude  Anomalies 

An  additional  measure  of  regional  attenuation  is  provided  by  amplitude 
anomalies  of  short-period  S  waves.  Short-period  S  waves  are  commonly  not 
observed  for  shallow  focus  events,  but  are  fairly  common  for  deep  earthquakes. 
This  fact  by  itself  puts  some  limitation  on  the  values  of  t*  that  we  shall 
discuss  later. 

If  regional  P  magnitude  anomaly  patterns  are  caused  by  variations  of  Q, 
then  regional  amplitude  anomaly  patterns  of  S  should  resemble  those  of  P. 

This  appears  to  be  the  case.  In  Figure  41  we  show  average  SH  trace  amplitude 
residual  terms  expressed  in  units  of  log^ (Amplitude)  from  the  seven  deep 
earthquakes  listed  in  Table  VI.  The  amplitudes  were  corrected  for  SH  radia¬ 
tion  pattern  using  the  double  couple  representation  of  the  sources.  Data 
points  close  to  nodal  lines  and  requiring  a  large  correction  were  omitted 
to  avoid  overcorrections  due  to  uncertainties  in  source  orientation.  The 
average  event  magnitude  was  used  for  normalizing  the  amplitudes  for  each 
event.  The  SH  residuals  clearly  show  patterns  similar  to  those  for  the 
P  waves  in  Figure  33  (Booth  et  al,  1974).  Furthermore,  as  Figure  42  shows, 
correction  for  crustal  amplification  effects  on  SH  using  flat  layered  models 
of  crust  under  LRSM  stations  does  not  remove  the  anomalies.  The  approach 
used  is  identical  to  that  presented  by  Der  et  al  (1979).  The  scatter 


-81- 


RKON 


M  tfl  cO 
O  rH 

O  <1)  r-t 

U  Cd 


0  73 
O  <U  73 
CUD 
Cd  <A  4-» 
P  cd 
o)  -nfl 

73  73  *H 
3  rt  u 
<M  (A 
•r^  73 
«H  C 
{X  cd  73 

i  CA  c3 
C 

d)  m  <i> 

>  a)  o 


Cd  4-> 

K  cx  w 

co  *H 

a  73 

73  O 
o  *«-<  * 

M  U  0) 

n  cd  o) 

<D  *H  73 
CU  73  P 
I  cd  u 
u  U  »H 

n  c 

O  U  toO 
&  O  cd 
w  w  0 


XJ  • 
4-1  c 
POO 
o  -W 

CO  C  U 

o  cd 
0J  «-l  P 
^  U  C 
U  3  D 
C  -M 

C  *d  ti 

•H  0  Cd 
•H 

»-i  73  U 
P  *H 

O  QJ  -u 
(0 

O  U  !fl 

(AMD 

o)  c  d 

p  cd 

r-i  4J 

cd  cd  >> 

>  U  ^3 

•w 

O  73  P 

>  c  o 

•H  *H  vl 

4J  00 

Cd  (A  CD 
W)  0)  P 
0)  J-> 

c  cd  co 

<u  U2  3 

00  4J 

• V 

HOC 
>-J  U  Tt 
•H 

C  09 
•  D  <U 
C  •« 

O  C  3 
M  W 
■u  4)  *H 

4)  4J  i— t 

u  co  a 
•w  m  e 
u  M 


83' 


is  still  quite  large,  a  fact  that  is  not  surprising  since  S  waves  are  also 
subjected  to  the  same  lateral  inhomogeneity  effects  as  P  waves,  but  the 
pattern  is  clear  nonetheless.  Unfortunately,  short-period  S  waves  are  not 
sufficiently  abundant  to  achieve  the  same  accuracy  in  the  estimation  of 
station  terms,  but  most  of  the  negative  residuals,  especially  the  large 
ones,  are  in  the  western  United  States. 

S  Wave  Spectral  Measurements 

According  to  most  evidence,  t*  "  4  t*  appears  to  be  valid  for  regional 
as  well  as  teleseismic  changes  in  t*  (Solomon  and  Toksoz,  1970;  Der  et  al, 
1980).  Consequently,  spectral  changes  in  S  should  be  quite  noticeable  from 
region  to  region.  This  appears  to  be  true  also.  In  Figures  43  to  49  we 
show  tracings  of  short-period  S  waves  at  various  LRSM  stations  across  the 
U.S.  Pointers  show  the  location  of  the  recording  LRSM  station  on  the  map, 
and  the  instrument  gains  are  given  above  each  trace.  Since  one  of  the  per¬ 
pendicular  horizontal  components  of  the  LRSM  stations  is  oriented  towards  NTS, 
the  tracings  were  done  on  the  component  closest  to  the  transverse  direction 
(SH)  to  the  event.  Although  a  few  of  these  components  may  be  misoriented  by 
as  much  as  45°,  the  figures  show  a  clear  tendency  for  the  stations  in  the 
southwestern  U.S.  to  have  lower  amplitudes  and  lower  frequencies.  These  are 
typical  examples  of  this  phenomenon  and  are  not  accidents  of  fault  plane 
orientation  and  directivity.  All  broad  band  S  waves  (sometimes  containing 
significant  signal  energy  up  to  2.5  Hz)  show  this  phenomenon  regardless  of 
fault  plane  orientation.  These  events,  not  previously  analyzed  in  our  recent 
study  (Der  et  al,  1980),  demonstrate  that  visible  regional  differences  in  the 
frequency  content  of  o  waves  are  the  rule  and  there  is  nothing  exceptional 
about  the  data  presented  by  Der  et  al  (1980).  We  have  inspected  S  waves  from 
large  numbers  of  deep  events  all  showing  similar  variation. 

Differential  attenuation  manifests  itself  in  various  ways  in  the  time 
domain.  Wide  band  signals  with  considerable  high  frequency  content  undergo 
a  significant  lengthening  of  dominant  periods  by  the  preferential  reduction 
of  the  high  frequency  end  of  the  spectrum.  Narrow  band  signals  containing 
less  high  frequency  show  little  change  in  waveform.  For  both  types  of  signals 
the  overall  trace  amplitude  is  reduced,  and  this  reduction  is  greater  for  the 
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Figure  43  Tracings  of  short-period  SH  phases  at  LRSM  stations 
across  the  United  States.  Depending  on  the  fre¬ 
quency  content  of  the  signals,  the  time  domain  mani¬ 
festations  of  anelastic  attenuation  vary,  but  the 
overwhelming  majority  of  the  signals  show  a 
diminution  of  amplitudes  and  the  decrease  of  high 
frequency  content  in  most  of  the  WUS  with  especially 
severe  effects  in  the  southwestern  United  States.  No 
corrections  for  radiation  patterns  were  made  in  these 
figures.  Instrument  gains  are  shown  on  each  trace. 
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44  Tracings  of  short-period  SH  phases  at  LRSM  stations 
across  the  United  States.  Depending  on  the  fre¬ 
quency  content  of  the  signals,  the  time  domain  mani¬ 
festations  of  anelastic  attenuation  vary,  but  the 
overwhelming  majority  of  the  signals  show  a 
diminution  of  amplitudes  and  the  decrease  of  high 
frequency  content  in  most  of  the  WUS  with  especially 
severe  effects  in  the  southwestern  United  States.  No 
corrections  for  radiation  patterns  were  made  in  these 
figures.  Instrument  gains  are  shown  on  each  trace. 
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Figure  46  Tracings  of  short-period  SH  phases  at  LRSM  stations 
across  the  United  States.  Depending  on  the  fre¬ 
quency  content  of  the  signals,  the  time  domain  mani¬ 
festations  of  anelastic  attenuation  vary,  but  the 
overwhelming  majority  of  the  signals  show  a 
diminution  of  amplitudes  and  the  decrease  of  high 
frequency  content  in  most  of  the  WUS  with  especially 
severe  effects  in  the  southwestern  United  States.  No 
corrections  for  radiation  patterns  were  made  in  these 
figures.  Instruments  gains  are  shown  on  each  trace. 
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Figure  47  Tracings  of  short-period  SH  phases  at  LRSM  stations 
across  the  United  States.  Depending  on  the  fre¬ 
quency  content  of  the  signals,  the  time  domain  mani¬ 
festations  of  anelastic  attenuation  vary,  but  the 
overwhelming  majority  of  the  signals  show  a 
diminution  of  amplitudes  and  the  decrease  of  high 
frequency  content  in  most  of  the  WUS  with  especially 
severe  effects  in  the  southwestern  United  States.  No 
corrections  for  radiation  patterns  were  made  in  these 
figures.  Instrument  gains  are  shown  on  each  trace. 
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Figure  48  Tracings  of  short-period  SH  phases  at  LRSM  stations 
across  the  United  States.  Depending  on  the  fre¬ 
quency  content  of  the  signals,  the  time  domain  mani¬ 
festations  of  anelastic  attenuation  vary,  but  the 
overwhelming  majority  of  the  signals  show  a 
diminution  of  amplitudes  and  the  decrease  of  high 
frequency  content  in  most  of  the  WUS  with  especially 
severe  effects  in  the  southwestern  United  States.  No 
corrections  for  radiation  patterns  were  made  in  these 
figures.  Instrument  gains  are  shown  on  each  trace. 
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Figure  49  Tracings  of  short-period  SH  phases  at  LRSM  stations 
across  the  United  States.  Depending  on  the  fre¬ 
quency  content  of  the  signals,  the  time  domain  mani¬ 
festations  of  anelastic  attenuation  vary,  but  the 
overwhelming  majority  of  the  signals  show  a 
diminution  of  amplitudes  and  the  decrease  of  high 
frequency  content  in  most  of  the  WUS  with  especially 
severe  effects  in  the  southwestern  United  States.  No 
corrections  for  radiation  patterns  were  made  in  these 
figures.  Instrument  gains  are  shown  on  each  trace. 
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wide  band  signals  (Der  and  McElfresh,  1980).  The  above-mentioned  tracings 
show  examples  of  both  types  of  behavior.  The  nature  of  signals  can  be  easily 
judged  "by  their  appearance  at  the  high  Q  shield  sites.  WWSSN  records  show 
similar  variations,  but  the  high  frequency  content  riding  over  lower  frequency 
waves  is  suppressed  somewhat  relative  to  those  shown  due  to  the  different 
instrument  response. 


The  differential  in  SH  wave  trace  amplitudes  reflects  the  regional  Q 
differential  to  a  lesser  extent  than  would  the  amplitude  at  the  dominant 
frequencies  of  the  signals  in  the  low  Q  regions.  Many  S  signals  in  the  north 
central  U.S.  have  dominant  frequencies  at  slightly  less  than  1  Hz;  such 
frequencies  appear  to  be  completely  absent  in  the  southeastern  U.S.  The 
Q  differential  can  effectively  be  measured  only  in  the  frequency  domain,  of 
course. 


In  a  previous  study  (Der  et  al,  1980)  we  have  shown  that  the  regional 

averages  of  t*  differentials  are  3  to  4  times  larger  than  the  average  t* 
s  P 

differentials. 

Outlining  the  Regional  Variations  of  Q  Under  the  United  States 

In  spite  of  the  many  uncertainties  with  regard  to  source  spectra,  source 

mechanism,  absolute  source  strengths  and  the  precise  determination  of  Q  in 

the  Earth,  reasonable  limits  on  the  values  of  Q  can  be  established  In  the 

short-period  band  by  utilizing  the  great  sensitivity  of  high  frequency 

energy  to  Q.  Consider  the  consequences  of  t*  ~  1  sec  and  t*  -  4  sec,  values 

P  8 

widely  used  in  long-period  simulations  and  also  claimed  to  be  valid  in  the 

short-period  band.  Figure  50  shows  spectra  of  P  and  S  waves  from  some  deep 

events  observed  in  the  north  central  United  States.  If  one  generously  allows 

for  the  source  depth  by  halving  the  t*  values  to  t*  *  0.5  and  t*  ~  2  sec  and 

0  -1  -2  s -3 

plots  the  expected  spectral  falloff  rates  of  u  ,  u>  ,  u>  ,  and  w  in  the 
source  spectra  (solid  lines),  there  is  a  sizable  discrepancy  at  the  high 
frequency  end  in  each  case.  This  discrepancy  is  especially  obvious  in  the  S 
waves  where  it  amounts  to  more  than  two  orders  of  magnitude.  Furthermore, 
for  S  waves  it  is  clearly  prevalent  even  in  the  vicinity  of  1  Hz.  The  ob¬ 


served  spectra  in  all  cases  are  incompatible  with  t*  ~  0.5  and  t* 
quire  considerably  lower  values  for  these  quantities. 


2  and  re- 
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Figure  50  Spectra  of  signal  and  noise  (upper  and  lower  dotted 

lines)  of  selected  F  and  S  waves  from  deep  earthquakes 
observed  in  the  north-central  U.  S.  The  spectra  were 
corrected  for  instrument  response.  The  falloff  rates 
of  theoretical  spectra  (solid  lines)  (assuming  various 
falloff  rates  for  the  source  spectrum  with  t*  •  .5  and 
t*  -  2.  and  also  allowing  for  the  source  deplh)  lead 
to  discrepancies  of  several  orders  of  magnitude.  These 
values  of  t*,  commonly  used  in  time  domain  work,  are 
therefore  unacceptable.  The  values  of  t*  ■  0.2  and 
t*  -  0.8  (dashed  lines)  with  an  at  fallSff  in  the  as¬ 
sumed  source  spectra  fit  better,  but  even  these  values 
may  be  too  high. 


We  estimate  from  a  larger  data  set  that  t*  <  0.2  and  t*  <  0.8  sec  is 

P  s 

needed  to  satisfy  these  data.  This  is  illustrated  by  the  dashed  lines  in 

-2 

Figure  50  assuming  a  falloff  rate  of  u  ,  which  is  the  most  likely  value  for 

large  events  such  as  these.  This  Implies  that  t*  £  0.2  and  t*  0.8  for 

P  8 

deep  events  and  that  for  shallow  events  in  this  type  of  upper  mantle  (shield)  at 

both  ends  of  the  paths,  t*  £  0.4  and  t*  <_  1.6.  This  result,  together  with  the 

long-period  result  of  t*  -  1  and  t*  ~  4  along  similar  paths  by  other  workers, 

p  s 

implies  that  t*  is  frequency  dependent.  We  must  point  out,  however,  that 
most  results  in  the  long-period  band  are  not  applicable  to  shields,  and  the 
possibility  remains  that  t*  is  low  but  not  dependent  on  frequency  for  shields. 


Consider  again  an  additional  constraint  that  limits  the  regional  varia¬ 
tion  of  t* — the  size  of  P  wave  magnitude  anomalies.  Since  these  are  of  the 
P 

order  of  0.3  magnitude  units,  the  corresponding  t*  variation  must  be  of  the 
order  of  0.2  sec  (Der  et  al,  1979).  This  is  the  same  result  that  one  arrives 
at  by  using  spectral  ratios  (Der  and  McElfresh,  1977;  this  report).  This 
implies  t*  **  0.6  for  paths  crossing  the  mantle  under  the  western  United 
States  and  terminating  in  a  shield  type  of  structure  at  the  other  end. 

Values  such  as  t*  =  1.3  (Hadley,  1979)  are  therefore  too  high  and  cannot  be 
accepted  in  the  short-period  band. 


As  far  as  measured  relative  regional  variations  of  short-period  t*  in 

3 

the  U.S.  are  concerned,  these  are  of  the  order  of  0.8  or  maybe  somewhat  less. 
The  existing  constraints  are  not  too  tight  on  these  (Der  et  al,  1980). 

t*  must  also  satisfy  another  condition  imposed  by  observations,  namely 
that  for  shallow  earthquakes  short-period  S  is  usually  not  seen.  This 
condition  is  quite  unspecific  since  it  is  not  clear  how  the  amplitude  of  S 
compares  to  that  of  P  in  the  short-period  band  at  the  source.  Assuming  the 
S  amplitude  to  be  five  times  that  of  the  P  (double  r  •«*-'  ),  one  gets  the 
result  that  for  t*  *  2  sec,S  will  have  5%  of  the  auq  >.  -  of  P  at  1  Hz.  A 

slightly  higher  t*  “  2.4,  corresponding  to  tectonic-to-shxeld  type  paths  and 
the  known  phenomenon  of  corner  frequency  shift  of  S  to  lower  frequencies, 
can  easily  account  for  the  observed  disappearance  of  S  into  the  P  coda. 
Therefore  our  values  appear  to  satisfy  this  condition  without  the  need  for 
t*  "  4,  a  high  value  commonly  claimed  by  some  researchers. 
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It  is  conceivable  that  one  could  occasionally  observe  a  t*  considerably 
larger  than  what  we  consider  representative  from  deep  events.  For  example, 
a  body  wave  could  encounter  pockets  of  attenuating  material  close  to  the 
source  (Sacks  and  Okada,  1974).  Such  observations,  however,  cannot  be  con¬ 
sidered  representative  for  the  mantle  under  the  north  central  United  States, 
a  common  medium  for  all  arrivals  observed  in  this  region.  Observations  of 
high  frequency  S  wave  energy  from  deep  events  would  be  unlikely  if  any 
observed  low  Q  is  attributed  to  the  mantle  under  the  observing  station, 
unless  we  presuppose  lateral  Q  variations  of  extremely  small  scale.  This 
is  because  the  raypaths  sample  a  relatively  small  region  under  the  station. 

All  of  the  data  presented  here  appear  to  be  compatible  with  an  apparent 

t*  of  0.1  to  0.2  second  for  shield  to  shield  type  paths  and  0.4  to  0.5  sec 
P 

along  paths  from  a  shield  to  the  WUS.  The  data  are  also  compatible  with 

the  measurements  for  explosions  listed  by  us  (Der  and  McElfresh,  1977)  and 

the  relative  t*  measurements  from  deep  events  (Der  et  al,  1980). 
s 

Correlation  with  Travel  Time  Delays  and  the  Extent  of  the  Mantle  Low  Velocity 
Layer 

It  is  not  necessary  for  travel  time  delays  to  correlate  with  regional 
variations  of  Q  since  seismic  wave  velocities  are  also  dependent  on  chemical 
composition.  Nevertheless,  the  correlation  between  the  regional  variations 
of  Q  outlined  above  and  travel  time  delays  reported  by  Sengupta  and  Julian 
(1976)  is  extremely  good.  The  major  features  of  their  results,  shown  in 
Figure  51,  include  large  delays  in  the  Basin  and  Range  and  SWUS  in  general, 
early  arrivals  in  the  shield,  and  slightly  late  arrivals  along  the  Atlantic 
seaboard  and  in  New  England.  The  late  P  arrivals  in  New  England  have  since 
been  studied  in  more  detail  (Taylor  and  Toksoz,  1979).  The  P  wave  travel 
times  at  0B2NV  are  1.7  sec  late  relative  to  RK0N  after  elevation  correction, 
and  because  of  this  NTS  fits  well  into  the  regional  pattern  presented  by 
Sengupta  and  Julian.  The  same  observation  can  be  made  with  regard  to  the 
S  wave  delays  in  North  America  as  compiled  recently  by  Wickens  and  Buchbinder 
(1980).  The  larger  S  delays  in  the  southwestern  corner  of  the  United  States 
appear  to  coincide  in  area  with  the  most  severe  attenuation  of  S  waves 
apparent  in  our  Figures  43  to  49.  There  are  also  indications  that  the 
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relationship  between  P  and  S  delays  is  different  in  the  EUS  versus  the  WUS 
(Romanowitz  and  Cara,  1980),  indicating  that  the  travel  time  delay  pattern 
cannot  be  associated  with  changes  in  thickness  of  simple  low  velocity  layers 
without  changes  in  velocity  contrast.  Low  velocities  in  the  upper  mantle 
also  manifest  themselves  in  the  diminution  of  phase  velocities  of  surface 
waves  with  significant  particle  motion  in  the  mantle.  There  appears  to  be 
a  correlation  between  low  surface  wave  Q  and  low  phase  velocities  worldwide 
(Lee  and  Solomon,  1979;  Nakanishi,  1979).  The  low  Q-low  mantle  velocity 
regions  are  also  characterized  by  negative  m^  residuals  (North.,  1977).  Low 
mantle  velocities  and  low  Q  thus  appear  to  be  areally  correlated  worldwide. 

Frequency  Dependence  and  Worldwide  Implications 

Strictly  speaking,  frequency  dependence  of  Q  can  only  be  established  if 
attenuation  is  measured  in  a  wide  frequency  band  using  short-  and  long-period 
instruments  over  the  same  paths.  Consideration  of  frequency  dependence 
introduces  a  large  number  of  additional  parameters  into  studies  of  attenuation. 
It  is  possible  that  the  form  of  frequency  dependence  changes  from  region  to 
region,  and  the  depth  distribution  of  Q  may  be  different  at  various  frequencies 
(Solomon,  1972;  Lundquist,  1979).  These  uncertainties  can  be  resolved  only 
by  further  detailed  studies.  A  large  number  of  reported  t*  determinations 
within  the  U.S.  in  the  long-period  band  were  obtained  by  time  domain  methods, 
but  these  cannot  be  accepted  as  valid  until  the  controversies  surrounding 
them  are  resolved.  We  have  shown  above  that  in  the  short-period  band  values 
of  t*  «  1,  as  often  obtained  by  time  domain  techniques,  are  unacceptable. 

In  the  following  we  shall  outline  possible  frequency  dependence  of  t* 
using  three  assumptions.  The  first  assumption  is  that  the  Q  and  velocity 
structures  are  basically  identical  for  all  shields,  so  that  absolute  attenu¬ 
ation  measurements  in  other  shield  regions  can  be  applied  to  the  north 
central  U.S.  The  second  assumption  is  that  t*  ”  4t*.  TM®  assumption  is 
supported  by  most  of  the  data  in  the  literature,  and  the  slight  modification 
due  to  losses  in  compression  proposed  by  Sailor  and  Dziewonski  (1978)  does 
not  affect  the  following  argument.  The  third  assumption  is  that  most  of 
the  losses  in  the  mantle  occur  in  the  upper  200  km  and  that  the  contribution 
of  any  low  Q  region  at  the  core  mantle  boundary  is  relatively  small.  All  of 
these  assumptions  agree  well  with  the  findings  of  research  to  date. 


-98- 


We  now  proceed  to  describe  attenuation  in  terms  of  t*  as  a  function  of 

P 

frequency  along  two  types  of  teleseismic  paths.  The  first  type  of  path 
crosses  a  shield-stable  platform  type  structure  at  both  the  downgoing  and 
upgoing  legs  of  the  path.  We  call  these  paths  "shield  paths"  in  the  following 
discussion.  The  second  type  crosses  a  WUS  (tectonic)  type  of  mantle  on 
either  leg  of  the  path.  We  term  these  paths  "shield-to-tectonic".  In  the 
previous  section  we  put  several  bounds  on  the  possible  values  of  t*  and  t* 
in  the  short-period  band  within  the  U.S.  These  are: 

.  t*  is  of  the  order  of  0.1  to  0.2  sec  for  shield  paths  in  the 
P 

0.5  to  4  Hz  range  (from  spectral  measurements). 

.  t*  ”  2  for  shield  paths  around  1  Hz. 

.  The  t*  differential  between  shield  and  shield-to-tectonic 
P 

paths  is  of  the  order  of  0.2  sec. 

These  constraints  are  supported  by  short-period  t*  measurements  from  other 

P 

regions.  Spectral  measurements  of  short-period  P  waves  having  shield  type 
paths  from  Asia  to  NORSAR  imply  t*  “  0.1  (Noponen,  1975,  Ringdal,  1976; 

Filson  and  Frasier,  1972).  The  spectral  differences  between  WUS  and  central 
Asia  earthquakes  and  explosions  at  NORSAR  imply  a  t*  differential  of  the 
order  of  0.2.  Spectra  of  teleseismic  P  waves  for  a  wide  variety  of  paths 
observed  on  shields  contain  significant  high-frequency  energy  in  the  3  to 
4  Hz  range,  essentially  ruling  out  any  constant  t*  ”  1  for  most  such  paths. 

Having  put  limits  on  t*  and  t*  within  the  short-period  band,  we  can 

proceed  now  to  review  the  evidence  in  the  long-period  band.  The  studies  by 

Solomon  and  Toksoz  (1970)  and  Solomon  (1972)  give  a  regional  t*  differential 

of  At*  ”0.5  sec  or  more  in  the  two  types  of  paths.  The  studies  of  Lee  and 
P 

Solomon  (1975  and  1979)  result  in  Qc  structures  that  imply  a  long-period  t* 

B  P 

differential  of  only  0.25.  Ray  tracing  through  the  Q  models  given  would 
yield  t*  of  the  order  of  0.6  to  1.7  in  the  eastern  U.S.  (EUS)  and  close  to 
unity  for  an  EUS-WUS  path.  This  by  itself  would  imply  frequency  dependence 
of  Q,  but  the  absolute  Q  values  in  these  models  are  rather  uncertain  due  to 
the  inherent  difficulties  of  measuring  Q  of  surface  waves  over  short  paths. 

In  any  case,  these  studies  indicate  that  the  upper  mantle  Q  also  varies 
regionally  in  the  long-period  band.  Therefore,  the  fitting  of  absorption 
band  models  that  do  not  allow  for  this  by  shifting  the  high-frequency  limit 
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only  (Lay  and  Helmberger,  1980)  unnecessarily  constrains  the  results.  The 
work  of  Nakanishi  (1979)  provides  further  indications  that  the  upper  mantle 
Q  measured  in  the  150  to  300  sec  period  range  is  high  under  shields.  It 
appears  from  his  work  that,  on  the  average,  anelastic  losses  under  shields 
are  less  than  those  associated  with  model  MM8  of  Anderson  et  al  (1965).  At 
teleseismic  distances  model  MM8  gives  a  t*  of  the  order  of  0.6  to  0.8  sec; 
thus  those  values  should  be  considered  as  upper  limits  of  t*  for  a  long- 
period  band  also. 


The  ideal  measurement  of  Q  under  shields  would  be  provided  by  multiple 

6 

ScS  phases.  Unfortunately,  there  are  no  studies  of  ScS  that  could  be  clearly 

associated  with  purely  shield  type  paths  (and  none  under  the  eastern  United 

States).  Nevertheless,  the  average  Q  ScS  values  for  the  whole  mantle  of 

B 

600  by  Kovach  and  Anderson  (1964)  and  580  by  Sato  and  Espinosa  (1967)  may  be 

indicative  of  Q  values  in  regions  above  the  downgoing  slab  in  South  America 

that  may  have  Q  characteristics  similar  to  shields  (Sacks  and  Okada,  1974). 

These  Q  values  are  considerably  higher  than  those  obtained  from  multiple 

ScS  studies  elsewhere  (Sipkin  and  Jordan,  1980),  but  the  corresponding  t*  at 

P 

teleseismic  distances  in  such  structures  would  still  be  0.4  to  0.5  sec,  twice 
the  apparent  t*  from  spectral  ratios  in  the  short-period  band.  Thus  even 
these  high  Q  values  imply  some  weak  frequency  dependence  of  Q  for  shield  type 
of  paths.  If  the  average  Q  under  shields  turns  out  to  be  lower,  as  suggested 
by  the  average  mantle  QScS  of  225  for  continents  (Sipkin  and  Jordan,  1980), 
the  frequency  dependence  would,  of  course,  be  stronger. 


The  above  constraints  allow  one  to  draw  the  preliminary  sketch  shown 

in  Figure  52  of  frequency  dependence  of  t*  for  the  two  types  of  paths  dis- 

P 

cussed  above.  The  most  natural  assumption  is  a  smooth  variation  that  appears 
to  be  supported  by  indications  of  smooth  changes  in  both  the  short-  and 
long-period  bands  (Archambeau  et  al,  1969;  Sato  and  Espinosa,  1967;  Brune, 
1977;  Yoshida  and  Tjusiura,  1975).  Such  a  gradual  change  does  not  greatly 
bias  any  relative  or  absolute  t*  measurements  from  spectral  ratios  although 
to  allow  for  such  bias  the  curves  are  drawn  higher  than  the  t*  determined 

from  spectral  ratios  assuming  a  constant  Q.  Details  of  the  regional  and 
frequency  dependence  of  Q  under  the  United  States  must  still  be  worked  out. 
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mixed  WUS— shield  path. 


We  ’-.ow  discuss  briefly  some  worldwide  measurements  of  Q  in  the  context 
of  frequency  dependence  and  regional  variations.  Worldwide  measurements  of 
Qgcs  311(1  ®antle  waves  (Sipkin  and  Jordan,  1979,  1980;  Nakanishi,  1979)  show 
that  the  regional  variation  of  Q  in  the  long-period  band  is  similar  to  that 
observed  in  the  short-period  band  (Bara2angi  et  al,  1975;  Oliver  and  Isacks, 
1967;  Solomon,  1972).  Attenuation  is  extremely  high  behind  island  arcs  in 
certain  areas,  under  mid-ocean  ridges  and  in  tectonic  regions.  It  is  low 
under  shields  and  old  ocean  basins.  This  appears  to  be  true  throughout  the 
0.003  to  4  Hz  band.  In  most  regions  of  the  Earth,  attenuation  measurements 

in  the  long-period  band  indicate  Q  values  that  would  give  t*  ”  1  and  t*  “  4. 

P  s 

Most  Q  models  derived  from  free  oscillation  data  also  imply  t*  ”  1  and 
t*  ~  4  (Anderson  and  Hart,  1977).  On  the  other  hand,  observations  of  high 
frequency  energy  in  the  3  to  5  Hz  range  are  quite  common  over  a  wide  variety 
of  teleseismic  paths  (Asada  and  Takano,  1963;  Takano,  1971;  Felix  et  al,  1971 
Noponen,  1975)  providing  the  observations  were  made  using  instrumentation 
peaked  at  high  frequencies  and  having  suitable  recording  systems.  A  large 
number  of  worldwide  estimates  are  given  by  Rivers  and  Der  (1980)  Indica¬ 
ting  t*  *  0.5  for  most  of  the  paths  studied.  Reports  of  similar  observations 
in  the  literature  are  too  numerous  to  quote  them  all.  There  is  thus  a  clear 
conflict  between  short-  and  long-period  Q  measurements  that  apparently  can 
be  resolved  only  by  assuming  a  frequency  dependent  Q  that  doubles  within  the 
range  of  0.01  to  2  Hz. 


CONCLUSIONS 


Large  scale  regional  amplitude  anomalies  exist  for  teleseismic  P  and  S 
waves  in  the  United  States.  The  regional  distribution  of  P  and  S  wave 
amplitude  anomalies  is  essentially  identical,  but  the  range  of  variation 
is  greater  for  S  waves.  The  amplitude  patterns  cannot  be  explained  by 
crustal  amplification,  and  corrections  for  crustal  amplification  leave 
a  pattern  of  high  amplitudes  in  the  north  central  United  States  and 
low  amplitudes  in  the  western  United  States — especially  in  the  south¬ 
western  part  of  the  country.  The  northeastern  United  States  is 
characterized  by  moderately  depressed  amplitude  levels. 

The  regional  amplitude  anomalies  correlate  with  spectral  changes  in 
both  P  and  S  waves;  low  amplitudes  are  accompanied  by  losses  in  high 
frequency  energy.  These  spectral  changes  are  quite  dramatic  in  S 
waves  from  deep  events  observed  across  the  United  States.  The 
regional  variation  of  teleseismic  t*  across  the  United  States  is  of 
the  order  of  0.2  seconds  for  P  and  about  three  to  four  times  that  for  S. 

The  existence  of  these  anomalies  confirms  the  hypothesis  that  the 
causes  of  these  variations  are  lateral  changes  in  Q  in  the  mantle  under 

p 

the  United  States. 

The  data  presented  indicate  that  mantle  attenuation  is  greatest  under 
the  southwestern  United  States,  including  the  Basin  and  Range  province, 
and  it  is  the  least  in  the  shield  region  of  the  north  central  Unitea 
States.  As  a  whole,  the  western  United  States  mantle  is  more  attenuating 
than  the  mantle  under  the  eastern  United  States.  The  northeastern  United 
States  appears  to  be  characterized  by  mantle  attenuation  greater  than  that 
of  the  shield  region  but  less  than  that  of  the  Basin  and  Range  province. 

The  SDCS  results  fit  well  into  the  regional  pattern  outlined. 

The  regional  variations  in  0  correlate  well  with  travel-time  residuals 
and  variations  in  the  upper  mantle  LVZ  that  are  derived  from  surface 
wave  studies  that  indicate  that  the  LVZ  is  also  a  low  Q  region. 

The  amount  of  high  frequency  energy  in  short-period  teleseismic  P  and  S 
waves  in  the  U.S.  and  worldwide  is  incompatible  with  the  values  of  t*  -  1  sec 
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and  t*  “  4  sec  derived  from  long-period  attenuation  studies  and 
s 

commonly  used  in  synthetic  simulations.  This  appears  to  indicate 
that  Q  is-  frequency  dependent  and  doubles  in  value  somewhere  in  the 
range  of  0.01  to  2  Hz. 
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PART  3 


BASIC  QUESTIONS  RELATED  TO  THE  ANALYSIS  OF  SHORT-PERIOD  DATA  AND  THE 
MEASUREMENT  OF  ATTENUATION  IN  THE  0.5  TO  5  HZ  BAND 

Interpretation  of  short-period  data  and  the  determination  of  Q  from  such 
data  presents  special  problems  that  are  of  no  concern  in  the  long-period 
band.  In  the  short-period  band  the  scale  of  inhomogeneities  in  the  crust 
and  the  surface  topography  become  comparable  to  the  wavelength,  and  that 
introduces  complications  due  to  focusing  of  seismic  energy,  crustal  amplifi¬ 
cations  of  signals  and  scattering.  In  addition,  the  frequency  content  of 
signals  becomes  a  critical  and  sensitive  determinant  of  Q  at  high  frequencies. 
The  tools  seismologists  use  in  analyzing  signals  are  the  measurement  of  wave 
amplitudes  and  spectra  and  the  matching  of  waveforms.  In  the  following 
sections  we  shall  evaluate  the  relative  effects  that  various  factors  have  on 
these  signal  characteristics.  We  shall  discuss  the  effects  of  Q  on  ampli¬ 
tude  and  spectral  measurements  and  critically  evaluate  time  domain  methods  of 
waveform  matching.  We  shall  also  estimate  the  possible  effects  of  instru¬ 
ment  nonlinearity  on  Q  measurements.  As  dictated  by  the  diversity  of  sub¬ 
jects,  this  part  of  the  report  consists  of  separate  sections  discussing  the 
above  topics.  These  are  referred  to  when  appropriate  in  the  previous  parts 
of  this  report. 
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Section  A:  The  Effect  of  t*  on  the  Absolute  Level  of  Spectra  in  the 
Short-Period  Band 

For  a  constant  t*  the  effect  of  attenuation  is  described  by  the  formula 

A  “  exp(-Trft*>  (1) 

where  A  is  the  wave  amplitude,  f  is  frequency  and  t*  *  T/Q  .  T  is  travel 

a  v 

time  and  is  c^e  average  quality  factor.  Let  us  examine  this  function 

for  some  values  of  t*  and  t*  frequently  mentioned  in  the  literature.  Figure 

p  s 

53  shows  a  plot  of  equation  (1)  for  various  t*.  As  the  figure  shows, 
t*  =  1  implies  that  4  Hz  amplitudes  are  reduced  by  a  factor  of  more  than 
P4 

10  relative  to  those  at  1  Hz.  Thus  the  assumption  of  t*  *  1  rules  out 
the  observation  of  4  Hz  energy  in  P  waves  for  all  practical  recording  systems 
presently  in  use  since  lower  frequencies  would  saturate  the  system  before 

4  Hz  energy  would  be  observable. 

Now  let  us  consider  the  consequences  of  various  constant  t*  on  the  spectra  of 

5  waves,  as  shown  in  Figure  54.  The  factor  of  4  increase  of  t*  relative  to 

t*  essentially  shifts  the  frequency  axis  by  the  same  factor  to  the  range 

around  1  Hz.  This  effect  is  even  more  severe  than  the  value  t*  =  3  sec 

s 

claimed  for  deep  earthquakes  by  Burdick  (1978)  and,  in  spite  of  the  counter¬ 
balancing  effect  of  most  short-period  instrument  responses,  effectively 
rules  out  any  observation  of  1  to  2  Hz  energy  in  S  waves  from  deep  events. 

It  appears,  however,  that  4  Hz  energy  is  routinely  observable  from  P 
waves  even  at  low  Q  sites  such  as  0B2NV.  This  is  demonstrated  in  Figure  55, 
where  we  show  five  typical  spectra  from  this  site.  Furthermore,  1  to  2  Hz 
energy  is  often  observed  in  S  waves  from  deep  earthquakes  as  shown  in 
Figure  50  of  part  2  of  this  report.  Observational  evidence,  therefore, 
precludes  the  general  use  of  such  high  values  as  t*  ’  1  and  t*  "  4  since, 
as  we  have  shown  above,  such  values  would  depress  the  higher  frequencies  by 
many  orders  of  magnitude  and  thus  render  them  unobservable. 

In  view  of  the  fact  that  spectra  of  short-period  body  waves  are  extremely 
sensitive  to  even  small  variations  of  Q.  we  must  conclude  that  studies 
claiming  these  high  values  are  incorrect.  We  shall  show  in  section  C  that 
distortion  of  the  shape  of  the  spectra,  especially  the  falloff  rate  toward 
high  frequencies,  is  also  primarily  determined  by  t*. 
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Figure  53  Diminution  of  P  wave  amplitudes  as  a  function  of 
frequency  for  various  values  of  t*. 
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Figure  54  Diminution  of  S  wave  amplitudes  as  a  function  of 
frequency  for  various  values  of  t*. 
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Figure  55  P  wave  spectra  at  OB2NV  shoving  significant  signal 
energy  at  4  Hz.  (All  these  spectra  have  a  minimum 
of  3:1  ratio  of  signal  to  noise  power.) 


Section  B: 


Time  Domain  Manifestations  of  Varying  t*  and  Their  Biasing 
Effect  on  the  Computation  of  m,  ~ ~~~ 


The  first  order  visible  effects  of  t*  in  the  time  domain  are  changes  in 
the  amplitudes  and  dominant  periods  of  the  body  waves  affected.  These 
changes  can  be  seen  to  various  degrees,  depending  on  the  spectral  character¬ 
istics  of  the  seismic  sources.  Figure  56  illustrates  this.  This  figure 
shows  a  wide  band  (left)  and  a  narrow  band  (right)  signal  causally  filtered 
by  various  t*.  The  narrow  band  Signal,  typical  of  large  events  used  in  time 
domain  simulations,  is  quite  insensitive  to  t*.  Note  that  its  waveshape  and 
dominant  period  change  little.  The  wide  band  signal,  on  the  other  hand, 
shows  a  drastic  change  in  general  appearance  with  increasing  t*,  and  the 
dominant  period  changes  from  0.5  sec  to  1.8  sec.  The  maxi-mum  trace  amplitude 
changes  by  a  factor  of  ten  for  the  wide  band  signal  and  by  a  factor  of  five 
for  the  narrow  band  signal  over  the  entire  range  of  0.8  sec  in  t*.  Thus  the 
effect  of  t*  on  two  time  domain  representations  can  be  visually  quite 
different,  and  not  all  signals  will  show  a  major  change  in  dominant  period. 


When  a  change  of  period  does  occur  it  will  introduce  a  paradox  if  m^  is 
computed  from  such  signals  using  the  standard  formula 

“b  "  lo8io  +  B(A8) 

where  is  the  trace  amplitude,  m(T)  is  the  magnification  of  the  instru¬ 
ment  at  the  dominant  measured  period  T,  and  B  is  a  distance-dependent 
correction  factor.  For  example,  if  m^  is  computed  for  the  wideband  signals 

associated  with  t*  *  0  and  t*  ■  0,6  at  the  left  of  Figure  56,  one  obtains  the 
P  P 

result  that  is  larger  for  t*  =  0,6.  This  is  the  reverse  of  what  one 
would  expect  on  physical  grounds  for  high  frequency  signals.  The  cause  of 
this  is  the  period  dependent  instrument  (LRSM  in  this  case)  factor  m(T)  that 
overcorrects  due  to  the  implicit  assumption  that  the  time  domain  amplitude 
measurement  is  associated  with  a  single  frequency.  This  paradox  affects  the 
body  wave  magnitude  and  amplitude  results  whenever  the  standard  procedure 
is  adhered  to.  For  moderate  changes  of  t*  the  actual  reversal  of  m^  values 
does  not  occur,  but  for  wide  band  signals  the  m^  procedure  tends  to 
de-emphaslze  the  attenuation  effect.  These  considerations  caused  us  to 
abandon  the  measure 


m 
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used  in  our  previous  reports  and  to  substitute  trace  amplitudes  A  in  our 
discussions. 

That  this  effect  is  not  a  hypothetical  presupposition  but  actually  does 
occur  is  shown  in  Figure  18  of  Part  I  and  the  accompanying  discussion 
(page  42  of  this  report).  The  fact  that  the  effect  of  attenuation  on 
critically  depends  on  source  spectra  makes  a  search  for  more  meaningful 
spectral  measures  of  body  wave  energy  necessary. 

These  remarks  were  made  in  order  to  point  out  how  deceptive  purely  time 
domain  observations  and  simulations  can  be.  These  problems  do  not  exist, 
on  the  other  hand,  if  comparisons  are  made  in  the  spectral  domain.  The 
advent  of  high  quality  recording  stations  with  large  dynamic  ranges  makes 
purely  time  domain  comparisons  obsolete,  since  spectral  calculations  are 
more  sensitive  and  reliable  indicators  of  variations  of  t*  with  frequency 
and  geographical  region.  Time  domain  comparisons  fail  to  utilize  the  broad 
band  information  available  in  signals  recorded  with  systems  of  high 
dynamic  range. 
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Section  C;  Various  Effects  on  Body  Wave  Spectral  Shapes  (Excluding 
Attenuation) 

Talcing  ratios  of  observed  body  wave  spectra  at  various  sites  in  order  to 
measure  relative  attenuation  and  taking  ratios  of  observed  spectra  to  model 
source  spectra  in  order  to  obtain  values  of  t*  are  widely  used  methods  in 
seismology.  These  techniques  utilize  the  slope  of  the  spectra  by  fitting 
straight  lines  to  spectral  ratios  on  a  semilogarithmic  plot,  discarding  the 
absolute  amplitude,  the  high  variability  of  which  is  unrelated  to  anelastic 
attenuation  (see  Section  D  of  this  report).  In  this  section,  we  shall  list 
and  evaluate  various  effects  not  related  to  t*  that  might  seriously  bias  the 
results  of  such  studies.  We  do  not  claim  that  these  effects  on  rare 
occasions  cannot  be  important,  but  we  merely  wish  to  consider  their  relative 
importance  compared  to  t*.  In  Section  A,  we  demonstrated  the  drastic  effect 
of  t*  on  body  wave  spectra.  In  the  following,  we  shall  demonstrate  that  the 
effect  of  other  factors  is  much  smaller. 

The  most  prominent  crustal  perturbation  of  spectral  slopes  at  the 
source  is  the  effect  of  surface  reflections.  Compared  to  these,  the  internal 
reverberations  have  a  small  effect  (Fuchs,  1966).  For  earthquakes,  this 
includes  the  phases  pP  and  sP,  the  relative  amplitudes  of  which  depend  on 
the  surface  reflection  coefficients  and  the  orientation  of  the  source 
mechanisms.  For  explosions,  the  pP  and  any  possible  spall  phase  can  affect 
the  spectrum.  However,  pP  and  the  corresponding  clear  spectral  nulls  are 
mostly  absent  in  the  observed  seismograms,  indicating  that  either  the  effec¬ 
tive  surface  reflection  coefficient  is  small  (due  to  either  scattering  or 
nonplanar  surfaces)  or  that  multipathing  effects  obscure  the  surface 
reflections.  In  any  case,  the  P  wave  spectra  may  be  affected  by  the  surface 
reflections. 

It  must  be  pointed  out  that  synthetic  results  based  on  elastic  flat 
layer  models  are  often  demonstrably  not  valid.  Such  calculations  generally 
predict  a  large  free  surface  reflection.  In  contrast,  the  actual  data  for 
most  explosions  do  not  show  such  large  secondary  arrivals,  and  the  spectral 
minima  (nulls)  to  be  expected  are  weak  or  not  detectable  (Der  and  McElfresh, 
1976).  It  appears  that  the  effective  reflection  coefficient  of  the  free 
surface  is  much  less  than  unity  at  most  places.  The  physical  reasons  for 


this  are  that  the  free  surface  Is  not  flat  and  that  there  exist  inhomogene¬ 
ities  close  to  the  source  that  distort  the  waveforms.  The  reduction  of 
effective  surface  reflection  coefficient  means  that  the  nulls  in  the  spectra 
disappear,  and  the  effect  on  the  t*  measurements  becomes  negligible.  Fur¬ 
thermore,  in  the  few  cases  when  these  nulls  are  apparent  in  the  data,  it  is 
easy  to  correct  the  spectra  for  pP  interference  and  eliminate  this  factor. 

For  earthquakes,  pP  and  sP  both  occur  with  amplitudes  varying  relative  to 
direct  P.  In  such  situations  clear  spectral  nulls  usually  do  not  appear. 

The  variants  of  such  simulations  may  include  the  spall  phase  or  a  pP  with 
the  spectra  modified  by  near-surface  layers.  In  any  case,  none  of  these 
models  can' consistently  mimic  the  effect  of  anelastic  attenuation  that 
consistently  suppresses  the  high  frequency  end  of  the  spectrum. 

Another  possible  effect  on  spectra  and  waveforms  is  that  of  multipathing 
and  focusing.  Studies  of  seismic  arrivals  at  arrays  reveal  that  secondary 
arrivals  of  energy  delayed  in  time  and  (often)  deflected  in  the  direction  of 
arrival  are  present  in  most  teleseismic  body  waves  (Mack,  1969) .  Some  of 
these  can  be  deterministically  modeled  by  an  uneven  Moho  or  deep  structures 
in  the  mantle  (Berteussen  et  al,  1975;  Capon,  1974;  Capon  and  Berteussen, 
1974;  Christofferson,  1975;  Dahle,  1975;  Dahle  et  al,  1975;  Haddon  and 
Husebye,  1978;  Hadley,  1979;  Chang  and  von  Seggem,  1980  and  many  others)  or 
by  the  theory  of  waves  in  homogeneous  random  media.  Th*  effect  of  identical 
multipath  arrivals  with  random  amplitudes  at  random  times  is  to  introduce 
fluctuations  in  the  body  wave  spectra,  but  it  cannot  introduce  a  consistent 
decrease  of  amplitudes  with  frequency  similar  to  the  factor  exp(-irft*). 

To  assess  the  effect  of  such  random  variations  we  computed  spectral 
ratios  between  individual  sensor  pairs  belonging  to  various  subarrays  at 
NORSAR  for  ten  events.  Since  we  do  not  believe  that  t*  actually  changes 
across  the  array,  the  observed  fluctuations  in  the  slopes  of  the  spectral 
ratios  between  0.5  to  4  Hz  must  reflect  the  random  effect  due  to  multi¬ 
pathing.  Figure  57  shows  the  histogram  of  slopes  of  spectral  ratios 
expressed  in  terms  of  apparent  t*.  The  standard  deviation  of  this  population 
is  0.06  sec,  showing  that  the  spectral  ratios  are  quite  stable  and  the 
scatter  is  quite  small.  Using  the  empirical  formula  Am^  a  1.35  At*  (Der  et 
al,  1979),  this  would  translate  into  Am^  -  0.08.  The  actual  variation  of 
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Figure  57  Relative  t*  between  subarrays  at  NORSAR  for  ten  tele- 
seismic  events.  The  histogram  shows  that  the  scatter 
of  At*  compared  to  that  of  amplitudes  is  small 
(o  *  0.06  sec).  This  illustrates  the  relative 
stability  of  spectral  measurements. 


Am^  measured  directly  is  about  0.4  magnitude  units  at  NORSAR.  This  demon¬ 
strates  that  spectral  slopes  are  more  stable  than  signal  amplitudes.  Large 
variations  of  amplitudes  across  LASA  were  also  reported  by  Chang  and  von 
Seggern  (1980).  These  are  generally  distance  and  azimuth  dependent. 

The  receiver  crust  is  another  factor  to  consider  when  evaluating  studies 
that  involve  spectral  ratios.  The  effect  of  the  receiver  crust  manifests 
itself  in  the  increase  of  the  peak  amplitude  of  the  signal  and  the  prolonga¬ 
tion  of  wavetrains  (ringing)  at  sites  with  thick  low-velocity  sedimentary 
cover.  In  the  frequency  domain,  the  slopes  of  amplitude  responses  of  most 
crustal  models  are  flat  when  fitted  with  least  squares  lines  over  the  0.5  to 
4.0  Hz  range.  Over  shorter  frequency  ranges  these  slopes  will  be  greater. 

In  general  the  presence  of  the  receiver  crust  does  not  seriously  affect  the 
spectral  ratios. 

A  few  remarks  need  to  be  made  about  detailed  slopes  of  crustal  responses. 

In  computing  synthetic  seismograms,  specific  crustal  structures  are  routinely 
used  to  model  the  crustal  reverberations  in  detail  in  the  time  domain.  This, 
on  the  surface,  conveys  an  impression  of  precision.  Experience  with  short- 
period  spectral  work  indicates,  however,  that  the  details  of  such  calculations 
both  in  the  time  and  frequency  domains  cannot  be  trusted.  An  indication  of 
this  is  that  for  short-period  waves  Phinney's  (1964)  radial  to  vertical 
spectral  ratio  method  never  worked,  the  probable  reason  being  that  at  most 
places  the  near-surface  geology  is  horizontally  heterogeneous  giving  differ¬ 
ent  spectral  ratios  in  various  directions.  The  dominant  effect  of  the  crust 
is  associated  with  acoustic  impedence  changes  due  to  variations  in  the  near-sur¬ 
face  materials  (Der  et  al,  1979)  that  causes  considerable  changes  in  the 
amplitude  levels  while  leaving  the  gross  spectral  slopes  unaffected. 

Another  objection  often  voiced  against  spectral  calculations  is  that 
scattered  high  frequency  energy  in  the  coda  can  significantly  bias  the  t* 
measurements.  First  of  all,  it  must  be  pointed  out  that  while  scattering 
near  the  receiver  may  apparently  enhance  high  frequencies  by  some  mechanism 
such  as  P  conversion  to  Rayleigh  waves,  such  scattering  cannot  generate 
high  frequency  energy.  Any  such  signal  energy  observed  must  have  come  from 
the  source.  The  simplest  test  to  assess  the  relative  contribution  of  the 
P-wave  coda  to  the  high  frequency  part  of  the  spectra  is  band  pass  filtering 
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of  P-wave trains.  Figures  58  to  60  show  some  examples  of  this.  In  filtering 
these  traces,  a  set  of  causal  band  pass  filters  with  a  flat  response  in  the 
band  indicated  on  the  figures  and  a  24  dB/octave  falloff  outside  the  band 
were  used.  Algorithms  for  causal  digital  filtering  are  described  in  many 
textbooks  (Oppenheim  and  Schafer  1975  for  example).  In  most  of  the  cases 
shown,  the  envelopes  of  the  wavetrains  are  essentially  similar  in  all  bands. 
This  indicates  that  our  use  of  a  9-second  signal  window  is  not  only  representa¬ 
tive  of  the  signal  spectra,  but  it  also  has  the  beneficial  effect  that  the 
spectra  are  more  stable  than  those  computed  from  shorter  windows.  Any 
enhancement  of  high  frequency  energy  in  the  coda  is  not  comparable  in  effect 
to  that  of  even  small  changes  of  t*.  Thus,  the  claim  that  the  scattering 
effect  is  significant  has  no  basis  in  fact. 

Finally,  an  effect  that  can  significantly  alter  both  the  observed  rela¬ 
tive  amplitudes  and  the  spectra  is  directionality  of  earthquake  sources. 
Although  considerable  work  has  been  done  in  the  long-period  band  to  model 
such  effects,  and  fair  success  has  been  achieved  in  modeling  waveforms,  not 
all  of  the  problems  have  been  solved  to  date.  Modeling  of  unequal  P  and  S 
corner  frequencies  is  still  deficient  (Molnar  et  al,  1973;  Hanks,  1980),  and 
source  time  functions  are  over-simplified.  In  the  short-period  band,  simu¬ 
lation  of  waveforms  is  largely  unsuccessful  for  earthquakes  and  doubtful  for 
explosions.  It  is  likely  that  short-period  radiation  patterns  are  extremely 
complex  for  large  and  moderate  sized  earthquakes.  While  at  long-periods, 
with  wavelengths  comparable  to  the  fault  length,  an  earthquake  may  resemble 
a  double  couple  source  with  some  directionality  component  added,  this  is  much 
less  likely  to  be  true  in  the  short-period  band  where  the  details  of  source 
motion  in  space  and  time  play  a  more  significant  role.  We  conclude  that 
modeling  of  large  earthquakes  in  the  0.5  to  5  Hz  frequency  range  is  beyond 
the  state-of-the-art  since  no  one  has  done  it  successfully.  In  the  absence 
of  modeling  capability  one  has  to  rely  on  averaging  over  many  events  to 
reduce  the  source  directionality  effects  in  order  to  estimate  path  charac¬ 
teristics.  Fortunately,  all  seismic  regions  contain  enough  variability  in 
source  mechanisms  to  make  this  possible. 

We  tabulate  the  main  conclusions  of  this  section  in  Table  VII. 
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10  sec  of  the  signal.  Therefore  taking  spectra  of 
the  first  9  sec  of  P  does  not  introduce  a  significant 
bias  in  t*  relative  to  that  computed  from  shorter 


TABLE  VII 


Section  D:  Perturbing  Effects  Influencing  Body  Wave  Amplitudes  (Excluding 
Attenuation) 

Some  discussion  of  amplitudes  was  presented  in  the  previous  section  where 
perturbations  of  the  spectra  were  compared  to  those  of  the  signal  amplitudes 
as  a  result  of  various  extraneous  factors.  In  this  section,  we  shall  discuss 
these  factors  in  more  detail  and  at  the  same  time  critically  examine  the 
proposition  that  averaged  P-wave  amplitude  levels  at  any  location  can  be  used 
as  the  sole  measure  of  attenuation  in  the  mantle  underlying  that  location. 

We  have  examined  the  effect  of  varying  constant  t*  on  amplitude,  and  it 
appears  that  the  measured  amplitude  differential  due  to  t*  variations  is 
dependent  on  the  source  spectrum.  It  is  larger  for  high  frequency  signals 
and  smaller  for  low  frequency  signals.  This  in  itself  makes  it  difficult  to 
relate  t*  directly  to  amplitude  residuals,  although  it  can  be  done  either 
empirically  (Der  et  al,  1979)  or  by  synthetic  simulation. 

The  effect  of  the  source  crust  can  be  reliably  modeled  if  the  source  is 
known  to  be  linear  elastic  and  if  the  crust  is  a  laterally  homogeneous  flat 
layered  medium.  In  most  cases,  neither  criterion  is  satisfied.  The  effect 
of  an  unknown  effective  surface  reflection  coefficient  on  the  time  domain 
amplitudes  can  be  quite  severe,  and  small  variations  in  the  reflection  coeffi¬ 
cient  and  the  source  depth  can  cause  large  changes  in  amplitude.  This  can 
be  avoided  by  picking  amplitudes  prior  to  the  arrival  of  surface  reflections. 

The  receiver  crust  has  a  large  effect  on  amplitudes.  At  sites  located  on 
thick  unconsolidated  sediments  the  body  wave  amplitudes  can  be  several  times 
larger  than  at  sites  on  solid  rock,  and  this  amplification  is  more  severe  for 
S  waves.  The  regional  geology  under  the  average  seismic  station  is  usually 
not  known  well  enough  to  model  such  effects  with  sufficient  accuracy, 
although  crustal  corrections  computed  by  flat  layered  models  achieved  a 
statistically  significant  reduction  of  variance  in  average  amplitude  measure¬ 
ments  for  a  larger  set  of  stations  (Der  et  al,  1979).  This  does  not  mean, 
however,  that  corrections  computed  for  individual  stations  are  accurate.  This 
is  not  surprising  since,  even  if  the  site  is  carefully  surveyed  with  refrac¬ 
tion  profiles  or  borehole  measurements,  the  shear  velocity-depth  distribution 
is  usually  unknown. 


Let  us  consider  now  the  effect  of  randomness  of  the  media  around  the 
sources  and  receivers.  Lateral  horizontal  inhomogeneities  around  sources  appear 
to  be  the  cause  of  broad  regional  systematic  variations  of  amplitude  that  can 
be  associated  with,  the  source  locations.  The  existence  of  such  effects  has 
been  verified  by  several  studies  (Hadley,  1979;  Butler,  1979),  and  this 
implies  that  in  order  to  measure  t*  on  amplitudes  this  random  effect  must  be 
removed  by  averaging  measurements  from  many  source  regions. 

The  effects  of  inhomogeneities  close  to  the  receiver  are  also  well  docu¬ 
mented.  These  effects,  well  known  to  researchers  involved  with  analysis  of 
array  data,  manifest  themselves  as  large  azimuthally  dependent  relative 
amplitude  variations  between  even  closely  situated  sites.  The  relative 
amplitude  patterns  are  repeatable  for  groups  even  at  similar  azimuths  and 
distances,  and  even  the  intersite  modifications  of  waveforms  are  repeatable 
and  describable  in  terms  of  transfer  functions  (Filson  and  Frasier,  1972; 

Chang  and  von  Seggern,  1980;  Butler  and  Ruff,  1980;  Lay  et  al,  1979).  These 
phenomena  were  recognized  early  in  array  work  and  have  been  analyzed  exten¬ 
sively  using  various  models  (Berteussen  et  al,  1975;  Capon,  1974;  Capon  and 
Berteussen,  1974;  Christof ferson,  1975;  Dahle,  1975;  Dahle  et  al,  1975; 

Haddon  and  Husebye,  1978;  Chang  and  von  Seggern,  1980  and  many  others). 

The  existence  of  focusing  phenomena  is  also  consistent  with  the  picture 
of  the  crust  provided  by  the  COCORP  studies  (Schilt  et  al,  1978),  which 
clearly  demonstrate  the  widespread  occurrence  of  lateral  Inhomogeneities  in 
the  crust.  It  is  not  unusual  for  sensors  spaced  only  a  few  km  apart  to  show 
fairly  large  amplitude  ratios.  This  hardly  comes  as  a  surprise  since  the 
geology  is  complex  enough  at  most  locations  and  such  effects  due  to  focusing 
and  defocusing  of  seismic  waves  are  to  be  expected.  It  was  found  by  Chang 
and  von  Seggern  (1980)  that  at  NORSAR  such  anomalies  tend  to  average  out  to 
zero  if  amplitude  measurements  at  a  large  range  of  distances  and  azimuths 
are  taken.  However,  this  cannot  be  done  at  most  stations  that  operate  for 
limited  lengths  of  time  due  to  the  non-random  geographical  distribution  of 
sources.  This  means  that  such  effects  cannot  always  be  removed.  Besides, 
there  is  no  guarantee  that  the  averjge  of  such  site  anomalies  approaches 
zero  in  all  cases,  even  if  equal  weights  are  given  to  a  wide  range  of 
azimuths  and  distances. 
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The  expected  consequence  of  the  near  source  and  receiver  inhonogeneities 
is  that  magnitude  levels,  even  at  closely  spaced  stations,  may  e^diibit 
differenced  that  cannot  be  explained  and  that  P-wave  magnitude  anomalies  at 
individual  stations  cannot  be  used  by  themselves  to  determine  the  degree 
of  attenuation  in  the  mantle  under  any  location.  The  attached  Table  VIII 
summarizes  the  effect  of  various  factors  Influencing  trace  amplitudes.  Most 
of  these  are  of  the  first  order. 
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TABLE  VIII 


The  Effect  of  Various  Factors  on  Absolute  Signal  Amplitudes 


Factor 

Anelastic  Attenuation 

Crustal  Amplification 

Surface  Reflections 

Scattering 

Multipathing 

Source  Directionality 


Near  Source  or 
Receiver  Focusing 


Properties 


Relative  Importance 


Changes  signal  amplitudes  First  order,  but  less  for  low 

frequency  (f  < 1  Hz)  signals 


Amplifies  signals 

Change  amplitude 

Unknown 

Unknown 

Unknown 


Large,  comparable  in  size  to  that  of 
t*  variations  to  be  expected.  Can 
be  estimated  to  some  degree 

Can  be  estimated  but  time  domain 
estimation  is  nonunique 

Unknown 

Unknown,  can  be  large 

Few  demonstrated  examples  in  the 
short-period  band,  no  adequate 
methods  exist  to  estimate  its 
importance  in  the  short-period 
band.  Can  be  removed  by  averaging 
over  many  events 


Changes  signal  amplitudes  Larger  than  that  of  t*  variations 

in  question.  Great  obstacle  to 
estimation  of  yields.  Absolute 
level  of  change  cannot  be  estab¬ 
lished  at  most  locations 
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Section  E;  A  Critique  of  Time  Domain  Methods 

A  large  number  of  papers  in  recent  literature  employ  time  domain 
matching  of  synthetic  waveforms  with  observed  waveforms,  during  which 
attenuation  (t*)  is  estimated  in  a  deterministic  multiparametric  scheme. 

The  common  approach  begins  with  the  source  by  computing  some  fault  imdel  or 
an  RDP,  then  continues  the  computation  of  the  waveforms  through  determinis¬ 
tic  models  of  the  crust,  upper  mantle,  receiver  crust  and  instrument 
response.  Several  free  parameters,  including  t*,  are  adjusted  in  the 
process  to  obtain  the  "best"  visual  fit  in  the  time  domain.  Those  values 
thus  obtained  are  assumed  to  constitute  a  valid  description  of  the  source 
and  path  properties. 

This  method  achieved  fair-to-good  success  in  modeling  seismograms  in 
the  long-period  band;  consequently,  researchers  were  emboldened  to  apply  the 
technique  to  short-period  data.  Unfortunately,  the  methodology  as  applied 
thus  far  has  so  many  flaws  that  most  of  the  results  should  be  declared 
invalid  in  the  short-period  band  and  doubtful  in  the  long-period  band.  We 
summarize  some  general  criticisms  below  (more  specific  objections  are  given 
elsewhere  in  this  report): 

1.  The  quality  of  the  matching  of  two  time  domain  waveforms  in  the  short- 
period  band  primarily  depends  on  the  low  frequency  end  of  the  spectrum  for 
large  events,  thus  ignoring  the  high  frequencies.  Waveforms  for  such  events 
are  extremely  insensitive  to  t*  (Der  and  McElfresh,  1980)  and  cannot  reliably 
measure  variations  in  t*  of  the  order  of  a  few  tenths  of  a  second  (See  Figure  56) . 

2.  The  practitioners  of  this  method  create  the  overall  impression  that  the 
mantle,  crust  and  source  parameters  are  well  known  and  precisely  controlled 
in  their  simulations.  We  submit  that  this  is  not  the  case.  Short-period 
data  are  characterized  by  large  spatial  fluctuations  in  waveshapes,  ampli¬ 
tudes  and  spectra  brought  about  by  small  scale  inhomogeneities  in  the  earth. 

Simple  parameterization  cannot  adequately  describe  these  fluctuations. 

3.  The  high  variability  and  scatter  inherent  in  short-period  data  requires 
statistical  techniques  for  the  extraction  of  meaningful  information.  In 
particular,  analysis  of  variance  of  the  parameters  in  the  problem,  such  as 
t*  in  our  case,  should  be  considered.  In  much  of  the  synthetic  work  no 


statistical  evaluation  of  the  stability  of  results  is  given,  and  the  data 
sets  fitted  in  the  short-period  band  are  often  unreasonably  small. 

4.  There  is  a  tendency  to  extrapolate  successful  theoretical  models  to  the 
short-period  band  or  to  use  theoretically  defined  but  restrictive  models  for 
simulating  the  data.  The  validity  of  each  type  of  theoretical  model  should 
be  independently  demonstrated  using  carefully  chosen  data  sets.  Only 
models  that  have  clearly  been  validated  should  be  used  in  fitting  parameters. 
Unfortunately,  this  is  frequently  not  the  case.  For  example,  the  Haskell 
type  of  propagating  fault  model  that  cannot  adequately  simulate  the  fre¬ 
quently  observed  inequality  of  P  and  S  corner  frequencies  is  used  uncriti¬ 
cally  in  many  simulations  (Molnar  et  al,  1969;  Hanks,  1980).  Another 
example:  fitting  Minster's  absorption  band  model  to  data  (Lay  and  Helmberger, 
1980)  unnecessarily  restricts  the  class  of  obtainable  solutions. 

5.  As  far  as  attenuation  studies  are  concerned,  a  considerable  amount  of 
detail  in  time  domain  modeling  is  irrelevant.  Since  we  are  interested  pri¬ 
marily  in  the  spectral  content  of  the  source  versus  that  of  the  observed 
waveform,  details  of  time  domain  waveforms  depending  mostly  on  phase  pro¬ 
perties  of  spectra  are  of  no  interest.  Simple  limiting  arguments  with  regard 
to  spectra  are  sufficient  to  put  reasonable  bounds  on  attenuation. 

6.  Time  domain  methods  often  end  up  with  the  same  t*  '  1  and  t*  ’  4  in  the 

p  s 

long-period  band  for  travel  paths  involving  a  variety  of  upper  mantle  struc¬ 
tures.  This  conflicts  with  the  observed  regional  variations  of  surface  wave 
attenuation  (Nakanishi,  1979;  Mills,  1978;  Lee  and  Solomon,  1979).  It  also 
appears  that  if  one  considers  only  the  attenuation  results  obtained  by  not 
using  the  time  domain  methodology,  these  correlate  well  with  the  regional 
patterns  of  surface  wave  and  short-period  body  wave  attenuation  in  the 
studies  just  quoted.  It  seems,  therefore,  that  results  of  time  domain 
studies  with  regard  to  attenuation  do  not  make  sense  in  the  broad  geophysi¬ 
cal  context,  which  indicates  that  the  methodology  is  not  suitable  for 
measuring  attenuation. 

7.  In  the  short-period  band,  the  t*  obtained  by  time  domain  methods  is 
often  demonstrably  wrong  (pages  105  to  107 of  this  report).  Furthermore,  t* 
reported  by  various  authors  far  similar  paths  using  the  time  domain  methods 


differ  greatly.  For  teleseismic  paths  to  NTS,  Burdick,  and  He’  berger  (1979), 
Hadley  (1979),  and  Helmberger  (1973)  use  t*  values  of  1.0,  1.5  ^..d  0.7 
respectively.  Besides  the  fact  that  two  of  these  values  are  clearly  impos¬ 
sible  (pages  93  and  108  of  this  report),  the  range  of  variation  is  such  that, 
if  true,  they  would  imply  an  variation  of  about  0.67.  Since  the  regional 

magnitude  variations  we  are  attesting  to  explain  are  of  the  order  of  0.2  to 
0.3,  such  instability  in  the  results  cannot  be  tolerated,  and  time  domain 
methods  as  practiced  today  thus  have  no  place  in  magnitude-yield  studies. 

We  must  point  out  that  spectral  methods  yield  short-period  t*  that  are 
repeatable  and  quite  close  in  value,  as  reported  by  several  authors  quoted 
in  the  literature  and  in  this  report. 

8.  A  variant  of  time  domain  waveform  matching  is  the  use  of  synthetics  in 

both  the  short-  and  long-period  bands  (Burdick,  1978;  Hadley,  1979).  This 

is,  in  effect,  a  variation  of  the  spectral  ratio  method.  Although  this 

approach  has  a  potential  for  usefulness,  no  conclusive  results  have  been 

produced  thus  far  that  can  be  applied  to  the  short-period  band.  Burdick's 

results  for  short-period  S  are  outside  the  frequency  range  of  interest  (the 

dominant  frequency  of  these  short-period  S  waves  is  0.25  Hz).  Hadley's 

result  of  t*  “  1.3  is  not  unique.  Der  and  Blandford  (1979)  showed  that 
P 

plausible  modifications  to  more  reasonable  values  of  insufficiently  con¬ 
strained  parameters  in  Hadley's  simulations  can  reduce  t*  to  0.6.  It  appears 
again  that  the  time  domain  amplitudes  in  these  simulations  are  too  sensitive 
to  the  various  unconstrained  factors. 

Summarizing  this  section,  it  appears  that  although  time  domain  methods 
have  contributed  a  great  deal  to  our  knowledge  of  earth  structure  and,  to 
a  lesser  degree,  of  source  mechanisms,  they  are  not  suitable  for  studying 
attenuation  unless  radical  changes  are  made  in  the  methodology.  For  the 
reasons  listed  above  and  in  the  discussions  elsewhere  in  this  report,  we 
shall  disregard  all  results  with  respect  to  attenuation  produced  by  such 
methods. 
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Section  F;  The  Possibility  of  Generation  of  High  Frequencies  in  the  Recorded 
Signal  by  Instrument  Nonlinearity 

In  a  recent  paper  Sacks  (1980)  suggested  that  the  high  frequency  energy 
observed  In  short-period  seismograms  may  not  be  actual  signal  energy  but 
rather  an  artifact  caused  by  nonlinear  distortion  in  the  system  when  a  high 
amplitude  but  low  frequency  P  or  S  wave  excites  it.  This,  if  true  for  the 
LRSM,  SDCS,  or  LASA-NORSAR  systems,  would  be  quite  serious,  especially  with 
regard  to  attenuation  studies.  Although  it  seems  implausible  that  work  over 
the  course  of  20  years  has  overlooked  such  an  effect,  and  although  recording 
non-linearities  are  routinely  recognized  and  avoided,  we  have  re-examined 
this  question  and  concluded  that  with  such  instrumentation  there  is  no 
serious  possibility  of  such  effects. 

We  have  analyzed  some  constant  amplitude  harmonic  (constant  current, 
equivalent  to  constant  ground  acceleration)  calibration  signals  that  were 
routinely  run  for  the  LRSM  system  at  high  equivalent  ground  motion  amplitude 
levels.  These  high  levels  resulted  in  the  FM  analog  high-gain  recording 
system  being  close  to  the  nonlinear  (clipping)  levels.  The  selected  station 
was  MNNV  for  31  August  1963.  We  computed  amplitude  spectra  on  the  calibra¬ 
tions  as  well  as  on  the  preceding  noise  background.  Some  examples  of  these 
are  shown  in  Figures  61  to  64.  The  calibration  frequencies  range  from  0.33 
to  5.0  Hz.  These  calibrations  were  run  routinely  on  the  last  day  of  each 
month  at  every  LRSM  station,  and  the  resulting  system  response  was  plotted 
in  the  logs.  The  signal  spectra  all  show  clear  peaks  at  the  center  frequency 
of  the  calibrations,  and  the  calibration  spectra  quickly  descend  to  noise 
level.  Side  lobes  above  noise  level  are  visible  only  for  1.0  and  1.5  Hz. 
These  signals  are  less  than  a  factor  of  two  below  the  visible  FM  clipping 
levels.  In  fact,  the  1.0  and  1.5  Hz  spectra  had  to  be  computed  from  the 
radial  channels  because  the  vertical  channels  (for  these  frequencies  only) 
did  clip.  Inspection  of  the  film,  however,  shows  no  clipping.  This  proves 
that  the  clipping  is  in  the  FM  system  and  not  in  the  instrument.  Clipping 
of  this  sort  is  routinely  recognized,  and  the  standard  procedure  in  these 
cases  is  to  use  the  low  gain  recording  channels. 

Another  approach  to  the  question  of  non-linearity  is  to  observe  signals 
having  equal  levels  of  1  Hz  but  different  levels  of  5  Hz  energy.  Such 
observations  would  seem  to  rule  out  the  possibility  that  the  5  Hz  energy  was 
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Spectra  of  steady  state  calibration  signals  at  LRSM 
stations.  The  spectra  are  dominated  by  the  frequency 
of  the  input  signal  with  only  minor  contamination  by 
harmonics  presumably  generated  by  nonlinearity. 
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Figure  62  Spectra  of  steady  state  calibration  signals  at  LRSM 

stations.  The  spectra  are  dominated  by  the  frequency 
of  the  input  signal  with  only  minor  contamination  by 
harmonics  presumably  generated  by  nonlinearity. 
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Figure  63  Spectra  of  steady  state  calibration  signals  at  LRSM 

stations.  The  spectra  are  dominated  by  the  frequency 
of  the  input  signal  with  only  minor  contamination  by 
harmonics  presumably  generated  by  nonlinearity. 
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Figure  64  Spectra  of  steady  state  calibration  signals  at  LRSM 

stations.  The  spectra  are  dominated  by  the  frequency 
of  the  input  signal  with  only  minor  contamination  by 
harmonics  presumably  generated  by  nonlinearity. 


created  by  the  1  Hz  energy.  Figures  65  and  66  are  two  events  as  recorded  at 
the  center  element  of  the  C3  subarray  of  NORSAR.  These  signals  show  no  indi¬ 
cation  of  clipping  or  nonlinear  distortion.  Low  t*  values  are  based  on 
spectra  such  as  that  seen  in  Figure  65.  The  spectra  as  seen  through  the 
system  response  are  nearly  flat  so  that  a  small  non-linear  side  lobe  cannot 
account  for  the  high  frequency  energy.  Compare  Figures  65  and  66  which  have 
about  equal  1  Hz  levels  but  which  differ  at  5  Hz  by  a  factor  of  about  30  in 
amplitude.  These  differences  in  spectra  cannot  be  associated  with  non-line¬ 
arity. 

Additional  evidence  is  offered  by  close-in  observations  of  NTb'  explo¬ 
sions  with  the  same  instrumentation.  Figure  67  shows  recordings  at  KNUT  and 
MNNV,  at  distances  of  about  200  and  300  1cm  from  explosions  of  1100  kt  (BENHAM) 
down  to  perhaps  less  than  1  kt  (BUTE0).  Note  that  the  smaller  BUTEO  event 
has  the  higher  frequency  and  the  recorded  amplitude  for  BUTEO  is  equivalent 
to  a  large-magnitude  earthquake  teleseismically,  while  the  BENHAM  amplitude 
through  the  system  response  is  100  times  larger.  Figure  68  shows  the  raw 
unsmoothed  spectra.  For  BUTEO  the  spectrum  is  nearly  flat.  We  see  directly 
that  in  BENHAM  as  compared  to  BUTEO  the  high  frequencies  are  dramatically 
absent  above  3  Hz.  Also,  for  frequencies  below  about  7  Hz  the  amplitude 
spectra  are  down  only  about  30  dB  even  for  BENHAM.  In  Figure  69  we  show  the 
spectral  ratio  of  these  two  events.  Note  that  it  is  in  excellent  agreement 
with  cube  root  scaling  theory,  ranging  from  a  ratio  of  over  1000  at  low 
frequencies  to  a  ratio  of  about  10  at  5  Hz.  For  these  events  the  FM  systems 
do  not  clip  because  resistors  are  added  in  the  circuit  as  necessary  to  keep 
the  voltage  in  range.  This  does  show  that  the  basic  seismometer  system  must 
be  linear  to  much  higher  amplitude  levels  than  those  encountered  in  tele- 
seismic  practice. 

Many  of  our  arguments  are  also  based  on  short-period  S  waves,  such  as 
those  shown  in  Figures  43-49  of  the  first  section  of  this  report.  First  of 
all,  the  most  severe  clipping  could  not  cause  a  change  in  the  dominant 
periods  of  the  signals  unless  the  nonlinearity  is  clearly  seen  (by  reversing 
the  peaks).  Even  this  becomes  impossible  if  one  considers  the  instrument 
gains  in  these  figures,  if  the  S  amplitudes  are  equal  everywhere  as  some 
claim.  The  instruments  with  the  highest  gains,  thus  having  the  highest 
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Figure  65,60  Spectra  of  two  events  at  the  C3  subarray  of  NORSAR. 

For  the  same  amplitude  level  at  1  Hz,  the  level  at 
5  Hz  differs  considerably.  If  nonlinearity  were  the 
source  of  5  Hz  energy  for  the  various  events  these 
levels  should  be  the  same. 
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Spectra  of  two  events  at  the  C3  subarray  of  NORSAR. 
For  the  same  amplitude  level  at  1  Hz,  the  level  at 
5  Hz  differs  considerably.  If  nonlinearity  were  “he 
source  of  5  Hz  energy  for  the  various  events  these 
levels  should  be  the  same. 
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Figure  67  NTS  explosions  recorded  at  KNUT  and  MNNV. 


potential  to  generate  extraneous  high,  frequencies  by  nonlinearity,  should 
show  the  shortest  S  wave  periods  and  the  greatest  high  frequency  content. 
However,  in  reality  the  reverse  is  true.  Due  to  the  offsetting  inequalities 
of  S  wave  amplitudes  and  gains  the  signals  are  recorded  at  roughly  the  same 
recording  levels  on  the  magnetic  tape,  and  no  clipping  is  detectable. 

Since  the  mechanical  system  of  seismograms  is  highly  linear,  most  non¬ 
linearity  comes  from  the  FM  tape.  (As  noted,  the  calibration  signals  did 
not  clip  on  the  film  that  bypassed  the  FM  recording. )  Both  high  gain  (in 
the  SWUS)  and  low  gain  (in  the  EUS)  should  have  about  the  same  amount  of 
distortion,  yet  the  frequency  contents  are  visibly  unequal. 

Additional  data  is  available  from  the  tests  performed  in  the  manufac¬ 
turing  plant.  As  an  example,  we  choose  the  model  4681  seismograph  consisting 
of  the  "small  Benioff"  vertical  seismometer  and  a  phototube  amplifer.  This 
system  was  widely  used  in  the  LRSM  program.  Figure  70  from  TR-63-55 
illustrates  the  linearity  of  the  recorded  amplitude  versus  input  shake  table 
displacement  at  1  Hz.  The  small  deviations  seen  could  as  well  be  due  to 
inaccuracies  in  the  shake  table  driving  motors  as  to  the  seismometer.  This 
shake  table  is  still  in  use  in  Garland.  The  largest  signals  were  nearly 
10  cm  peak-to-peak  on  the  graphic  display,  and  no  other  frequencies  were 
visible  even  though  a  2  mm  ripple  could  have  been  seen.  This  places  the 
non-linearity  at  700  microns  more  than  34  dB  down.  700  microns  is,  of 
course,  much  larger  than  any  teleseismic  signal.  Figure  71  shows  a  complete 
frequency  response  for  this  instrument  on  a  shake  table.  Between  the 
plotted  points  there  is  not  much  room  for  a  spurious  resonance.  Other  raw 
data  plots  of  phase  and  amplitude  response  from  shake  tables  in  TR-59-14 
have  about  two  times  as  many  points  and  still  show  no  resonances.  The 
responses  are  generally  taken  at  an  input  displacement  of  9.2  microns  which 
is  much  larger  at  these  frequencies  than  most  teleseismic  signals. 

We  have  also  examined  the  develocorder  recordings  of  single  TFO  elements 
of  an  03  September  1967  event  off  the  coast  of  Peru.  This  is  shown  in 
Figure  72.  The  large  LR  waves  with  periods  of  14  seconds  as  seen  on  the 
short-period  instrument  arrive  around  21:42  and  are  on  scale,  with  an 
amplitude  of  15  cm  peak-to-peak  on  the  viewer.  There  is  also  visible  3  to 
5  Hz  energy  from  local  sources  of  noise  with  amplitudes  of  1  mm  peak-to-peak. 


-140- 


This  same  energy  is  visible  before  the  arrival  of  P.  In  this  case,  there  is 
no  apparent  non-linearity  of  the  seismograph  system  to  a  level  of  44  dB.  In 
terms  of  non-linearity  with  respect  to  ground  displacement,  taking  account 
of  the  differing  magnifications  at  14  seconds  and  3  Hz,  the  linearity  is  on 
the  order  of  44  +  80  =  124  dB.  In  this  case,  the  LR  also  includes  a  non¬ 
vertical  component  of  motion  that  eliminates  the  possibility  of  cross  talk 
between  components.  By  analogy,  the  short-period  (f  -  1  Hz)  oscillations  in 
the  S  waves  shown  in  Figure  72  cannot  be  caused  by,  say,  a  huge  0.1  Hz  S 
wave  exciting  the  seismometer.  Such  a  signal  would  still  be  dominantly  0.1 
Hz  in  frequency  as  seen  through  the  short-period  instruments. 

Further  arguments  showing  that  nonlinearity  is  not  an  admissible  explan¬ 
ation  for  low  t*  and  that  high  frequencies  are  an  integral  part  of  a  real 
signal  can  be  summarized  as  follows: 

a)  The  high  frequency  content  of  short-period  signals  recorded  on  identical 
systems  varies  with  geological-geophysical  setting. 

b)  High  frequency  signals  show  other  meaningful  geophysical  detail  such  as 
surface  reflections  of  pP,  sP  and  travel  time  triplications. 

c)  The  high  frequency  content  of  P  and  S  waves  differs  for  comparable  long 
or  intermediate  frequency  amplitudes. 

d)  The  high  frequency  content  of  signals  in  identical  geophysical  station 
settings  is  of  the  same  order  of  magnitude  regardless  of  component,  make 
and  construction  of  the  instruments  and  systems,  including  hydrophones. 

In  summary,  while  we  do  not  disagree  with  Sacks  that  nonlinearity  can, 
and  actually  does,  introduce  harmonics  at  levels  far  below  that  of  the  pre¬ 
dominant  signal,  we  must  conclude  that  in  our  analyses  such  effects  do  not 
play  a  significant  role. 


IMPLICATIONS  OF  THE  FINDINGS  OF  THIS  REPORT  TO  YIELD  ESTIMATION 

There  is  some  confusion  about  the  interpretation  of  the  results  of  the 
SDCS  experiment  and  the  significance  of  the  "discrepancy"  between  the 
bias  estimates  derived  from  amplitude  measurements  and  spectral  differences. 
There  also  appears  to  be  a  strong  feeling,  shared  by  some,  that  the  m^  bias 
derived  from  amplitudes  is  the  true  measure  of  bias  due  to  attenuation  because 
of  "reciprocity".  We  feel  that  a  few  clarifying  remarks  are  necessary. 
Beginning  with  the  reciprocity  theorem  as  it  is  stated  in  textbooks,  we  can 
say  that,  if  we  have  a  dilatational  source  and  a  dilatation  measuring  de¬ 
vice,  then  in  any  anisotropic,  inhomogeneous  medium  with  arbitrary  geometry, 
we  can  interchange  the  source  and  the  receiver  and  obtain  exactly  the  same 
seismogram,  if  the  source  is  represented  by  the  same  force  system  and  has 
the  same  time  function.  This  is  clearly  not  a  practical  situation,  and  we 
shall  now  explore  by  more  practical  examples  how  this  reciprocity  breaks 
down  as  we  move  away  from  the  above  mentioned  situation.  Let  us  assume  two 
granite  bodies  in  two  regions  and  explode,  underground,  two  nuclear  devices 
of  the  same  yield  at  the  same  depth  (but  not  at  the  same  time).  We  then  re¬ 
cord  the  seismic  waves  at  the  surface  near  the  explosion  sites  with  a 
vertical  seismograph.  Although  the  strict  mathematical  reciprocity  has 
already  broken  down,  we  can  still  expect  to  see  very  similar  seismograms 
since,  assuming  that  the  source  media  are  similar,  the  source  time  functions 
are  the  same  and  the  near  source  and  receiver  focusing  effects  are  still 
approximately  reciprocal.  If  we  assume  now  that  one  of  the  sites  has  a 
different  source  medium,  we  cannot  expect  to  see  the  same  seismogram,  because 
nuclear  devices  of  the  same  yield  do  not  give  the  same  equivalent  elastic 
source  strength  and  time  function  in  different  source  media.  The  amplitudes 
and  the  waveforms  at  the  two  sites  will  not  be  the  same.  Now  if  the  events 
are  moved  away  from  the  recording  sites  by  only  a  few  tens  of  kilometers, 
the  near  source  focusing  effects  will  be  different,  and  even  if  the  devices 
are  exploded  in  the  same  media,  the  seismograms  can  be  quite  different  in 
amplitude.  (The  seismograms  may  not  differ  in  waveform,  since  near  source 
focusing  probably  does  not  affect  them  as  much.)  The  next  step  is  to  move 
the  sites  far  away  from  the  explosions.  Even  in  identical  source  media  and 
with  identical  yields,  there  can  be  quite  large  differences  in  waveforms  and 
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amplitudes,  because  both  near  receiver  and  source  focusing  effects  will  be 
quite  different.  If,  as  in  the  SDCS  experiment,  we  use  sources  of  uneven 
azimuthal  distribution,  as  we  were  forced  to  do  because  of  the  given  limi¬ 
tations  of  global  seismicity  to  determine  differences  in  level,  a  "dis- 
crepancy"  between  the  measured  t*  and  residuals  is  not  surprising  but  to 
be  expected,  because  near  receiver  focusing,  which  is  azimuthally  dependent, 
cannot  be  averaged  out  entirely  with  the  available  data.  (We  must  note  that 

certain  types  of  inhomogeneities  that  consistently  focus  or  defocus 
teleseismic  arrivals,  there  is  no  guarantee  that  even  azimuthal  weighting 
would  help.)  The  only  way  the  measured  m^  would  be  a  valid  measure  of  "bias" 
including  the  Q  and  focusing  effects  is  if  we  had  a  measured  source  located 
where  the  observing  stations  are  located.  Therefore,  the  claim  that  "recipro¬ 
city"  requires  that  the  bias  measured  be  the  true  measure  of  anelastic 
loses  is  fallacious,  because  stations  are  not  located  near  the  event  epi¬ 
centers. 

The  relative  m^  levels  and  their  confidence  limits  resulting  from  this 
experiment  mean  relatively  little  for  individual  stations.  They  are  probably 
biased  by  focusing  and,  if  we  had  events  with  more  even  azimuthal  distribution 
(which  would  have  required  many  years  to  accumulate),  the  results  could  have 
been  quite  different  and  the  mean  m^  residuals  could  even  be  outside  the 
confidence  limits  given  in  this  report.  Some  examples  of  these  anomalies, 
are  the  .15  m^  difference  between  HNME  and  IFME  (which  cannot  be  adequately 
explained  by  near  surface  geology),  and  the  differences  between  the  FANV-GBNM 
pair  and  the  rest  of  the  WUS  stations,  although  the  t*  are  similar  within 
both  groups.  Focusing  provides  a  likely  explanation  for  these  differences 
and  the  references  in  this  report  demonstrate  that  this  phenomenon  is  real 
and  widespread  and  affects  seismic  waves  in  all  areas  of  the  world.  Common 
sense  dictates,  therefore,  that  while  we  should  not  disregard  amplitude  data, 
we  should  interpret  it  with  caution,  and  rely  more  on  spectral  data  to 
determine  anelastic  losses.  Anomalies  of  m^  are  more  meaningful  if  regional 
averages  which  can  be  interpreted  in  terms  of  Q  are  taken  (Der  et  al.,  1979; 

Lay  and  Helmberger,  1980). 

After  establishing  ,t1  it  the  SPCS  does  not  mean,  let  us  discuss  what 
was  actually  measured.  Since  spectra  are  less  affected  by  focusing  and  other 
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site  effects,  we  claim  that  the  contribution  of  Q  to  the  decrease  of  signal 
amplitudes  from  NTS  was  measured  by  estimating  differences  in  t*  relative  to 
a  shield  area.  Although,  according  to  the  best  evidence,  the  Soviet  test 
sites  appear  to  be  on  shields  or  stable  platforms,  the  findings  of  this  reporc 
do  not  bear  on  whether  RKON  or  HNME  is  a  better  analog  to  Kazakh  or  Novaya 
Zemlya.  Such  decisions  should  be  based  on  facts  involving  seismic  waves 
originating  from  those  areas.  At  the  SDAC  and  at  large  arrays,  many  years 
of  work  on  short-period  body  wave  spectra  support  the  idea  that  the  mantle 
under  shield  areas  is  less  attenuating  than  the  mantle  under  tectonic  regions. 
Nevertheless,  considerably  more  work  is  needed  to  test  the  validity,  with 
respect  to  mantle  attenuation,  of  possible  analogs  around  the  world. 

The  measured  m^  residuals  in  the  SDCS  experiment  are  in  rough  agreement 
with  the  t*,  but  no  exact  agreement  with  the  available  data  can  be  expected, 
since  even  the  time  domain  features  agree  with  the  interpretation  in  terms 
of  Q.  We  assume  that  the  available  data  was  sufficient  to  average  out  most, 
but  not  all  of  the  focusing  effects.  This  is  reassuring,  since  it  would  be 
quite  disturbing  if  the  relative  signs  of  the  m^  and  t*  residuals  were  not 
consistent  with  the  Q  interpretation.  Additional  S  wave  information  and 
broad  regional  studies  mentioned  in  this  report  also  confirm  that  the  Q 
effects  are  indeed  real. 

Although  body  wave  amplitude  measurements  by  themselves  are  not  suitable 
to  measure  attenuation  for  yield  estimation,  one  must  ultimately  use  wave 
amplitudes  for  estimating  yields.  First  of  all,  one  should  make  distinction 
between  the  t*  effect  and  the  effects  of  crustal  structure,  both  near  re¬ 
ceiver  and  near  source  focusing.  The  m^  bias  due  to  attenuation  alone  in 
upper  mantle  under  NTS  appears  to  be  in  the  0.20  to  0.2"  magnitude  unit 
range,  if  one  uses  the  empirical  formula  Am^  ~  1.35  At*  for  the  NTS  stations. 
This  slightly  higher  than  the  actually  measured  m^  differentials.  The 
formula  can  also  be  justified  on  the  basis  of  synthetic  studies  of  pulses 
from  nuclear  explosions.  Figure  73  shows  that  a  line  with  a  slope  of  1.35 
fits  curves  derived  from  synthetics  of  m^  versus  t*  quite  well.  This  figure 
also  shows  that  ir.  the  10  kt  to  100  kt  range  some  of  the  effects  on  m^  due 
to  varying  periods  discussed  in  Section  B  are  not  important,  although  they 
did  affect  the  earthquakes  used  for  measuring  Am^.  For  estimating  yields 
throughout  the  world,  upper  mantle  Q  should  be  thoroughly  mapped  under  the 
source  regions  of  interest  and  under  the  stations  used  for  determining  the 
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yields.  Crustal  effects  under  these  stations  should  also  be  estimated  and 
corrections  for  these  effects  should  be  made. 

Effects  due  to  near-source  focusing  can  be  determined  only  for  very 
small  source  regions  for  explosions  with  known  yields.  In  the  absence  of 
such  information,  averages  of  m^  for  a  large  range  of  epicentral  distances 
and  azimuths  will  probably  eliminate  such  effects.  In  general,  near-receiver 
focusing  cannot  be  eliminated  unless  all  receiving  stations  have  many 
sensors,  thus  enabling  one  to  outline  the  causative  structures  and  to  derive 
deterministic  formulas  for  corrections.  Since  the  requirement  for  the 
elimination  of  near-source  focusing  also  demands  a  large  range  of  azimuths 
and  distances  for  stations  around  the  source,  this  calls  for  many  large 
arrays.  We  do  not  consider  this  a  very  practical  or  economical  alternative, 
and  one  is  again  reduced  to  using  network  averages  of  m,  . 
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List  of  Events  Used  in  the  SDCS  Project  Along  With 
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APPENDIX  B 


Histograms  of  Magnitude  Differentials  for  Various 
Pairs  of  SDCS  Stations 
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Histograms 
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APPENDIX  D 


Histograms  of  Differentials  At*  for  Various  Pairs  of 
SDCS  Stations 
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APPENDIX  E 

Histograms  of  Dominant  Period  Differentials  AT  for 
Various  Pairs  of  SDCS  Stations 
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MEMORANDUM 


FROM: 


DATE: 


All  Recipients  of  VSC-TR-81-14 
Z.  A.  Der  and  T.  W.  McElfresh 
26  October  1982 


SUBJECT:  Explanation  of  amplitude  data  in  VSC-TR-81-14 

The  raw  station  data  presented  in  Appendix  A  of  the  final  SDCS 
report  (VSC-TR-81-14)  has  the  format: 


IF-ME  77.2  11.0  0.9 


Arrival 
H  M  S 

11  45  6.7 


where 


station  name 


station  to  event  distance  (deg) 

maximum  signal  amplitude*  (first  three  cycles) 

dominant  period  of  corresponding  cycle  (sec) 


minute 

second 


Arrival  Time 


gain  factor  (see  below) 
b-phase  amplitude,  if  measured 
b-phase  period,  if  measured 


*Attention  is  called  to  the  Amplitude  and  GAIN  columns.  If  no  gain 
factor  is  given,  then  the  amplitude  A  is  understood  to  be  in  units  of 
nanometers .  This  is  the  case  for  those  stations  having  digital 
recording  equipment . 

If  a  gain  factor  is  listed  for  a  station,  then  the  amplitude  A  is  in 
millimeters  as  messured  directly  from  the  film  viewer.  In  this  case, 
the  amplitude  in  nanometers  can  be  computed  using  the  equstion: 

A  -  500  A  /GAIN 
nm  mm 

This  conversion  was  applied  by  the  computer  program  that  produced  the 
histograms  of  magnitude  differentials  Am.  snd  trace  amplitude 
differentials  AAtr  shown  in  Appendices  B°and  C,  Thus  the  histograms 
and  results  of  thereport  are  valid  as  stated . 


The  few  repeated  events  in  Appendix  A  are  due  to  the  arrival  and 
processing  of  various  subsets  of  the  data  at  different  times.  This  also 
accounts  for  the  events  not  being  listed  chronologically. 

The  authors  regret  any  confusion  that  this  may  have  caused  our 
readers.  We  are  planning  to  re-issue  the  data  shown  in  Appendix  A  in  a 
unified  format  in  the  near  future. 
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